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Q I We study the possible holographic connection between the Regge limit in QCD and critical 
gravitational collapse of a perfect fluid in higher dimensions. We begin by analyzing the 
^ ' problem of critical gravitational collapse of a perfect fluid in any number of dimensions and 
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holographic correspondence between this process and the phenomenon of parton saturation 
in high-energy scattering in QCD is proposed. This holographic connection relates strong 
gravitational physics in the bulk with (nonsupersymmetric) QCD at weak coupling in four 
dimensions. 
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1 Introduction 

Recently it has been proposed that critical black hole formation might provide the gravitational 
dual description of diffractive scattering in QCD p[] (see also [21 [3l H] for other attempts at a 
holographic description of Pomeron physics). More precisely, the proposal identifies the onset 
of strong curvature effects in the five-dimensional scalar field gravitational collapse with the 
corresponding phenomenon in the BFKL evolution, leading to unitarity saturation in QCD 
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in the Regge limit. The computation of the corresponding exponents on the gravity and the 
QCD side render numerical values that agree quite well. 

There are still a number of open questions to be clarified in this conjectured duality. One 
of them is the reason to consider a scalar field instead of any other kind of collapsing matter. 
On physical grounds it seems that a perfect fluid would be more appropriate to describe 
the gravitational dual of hadron scattering. In the center of mass, in the Regge limit, the 
colliding protons are composed of a large collection of partons, in particular gluons. Hence, 
a hydrodynamic description, specially when we are interested in collective properties of the 
scattering process seems appropriate. In this respect it is important to notice that in four- 
dimensional gravitational collapse one finds that the critical exponents for a collapsing massless 
scalar field, and for a radiation perfect fluid are numerically very close. The best current value 
for the Choptuik exponent for a massless scalar field in four dimensions is 7 = 0.374 [5J, while 
for the conformal perfect fluid one has 7 = 0.3558 [HlllOj]. In five dimensions, the best value for 
a massless scalar is 7 = 0.431 ± 0.0001 [5J, while in this paper we compute that for a perfect 
conformal fluid in d = 5 the corresponding value is 7 = 0.388. The value to be compared with 
in the Regge limit of four- dimensional QCD according to [Ij is 7 = 0.409. All these numbers 
are reasonably "close" , but not close enough that further analysis are not required to establish 
that they are indeed related by some underlying physics. This implies a double task. On one 
hand it is necessary to draw a clear enough holographic map between the gravity and the gauge 
theory descriptions. On the other hand a more precise understanding of the QCD exponents 
is required. 

This paper is divided in two clearly distinct parts. The first part studies a well-posed 
problem: the determination of the critical exponents in the gravitational collapse of a perfect 
fluid in dimensions other than four, in particular in five dimensions. The second is dedicated to 
outline the possible holographic map of this critical collapse with the Regge region. Needless 
to say, this latter part is wanting a more concrete and straightforward formulation. 

The discovery of critical black hole formation has been one of the most exciting develop- 
ments in numerical general relativity (see |6] for a review). This phenomenon was discovered 
by Choptuik [7] when studying numerically the spherically symmetric collapse of a massless 
scalar field. Considering a set of initial conditions parametrized by a real number p, such that 
for "large" p a black hole is formed, whereas for "small" p the initial scalar field configuration 
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disperses to infinity, it was found numerically that there is a critical value p* marking the 
threshold of black hole formation. Surprisingly, for p ^ p* the size of the black hole follows a 
scaling relation, 

tbh ~ (p - P*)'^- (1-1) 

The critical exponent 7 is independent of the particular family of initial conditions chosen. It 
changes, however, when the scalar field is replaced by other types of matter. 

Notice that the critical solution corresponds to the formation of a "zero mass" black hole. 
If the scaling behavior (11. ip were observed only in a very restricted value of the parameter p 
close to p*, the very scaling law would be quite questionable. One could argue that higher 
curvature corrections to the Einstein equations might wash out this behavior. Indeed, the 
curvature at the black hole horizon behaves roughly like 1/m^ hence, as the mass goes to zero, 
the horizon size vanishes and the curvature grows quadratically. Therefore we should expect 
not only higher curvature effects, but also quantum gravitational effects to become relevant. 
However, the beauty of Choptuik's result lies in that the scaling law is observed in a long 
range of the parameter p before the black hole forms. It sets in before the curvature becomes 
overwhelmingly big. This means that we can expect that there is a region in parameter space 
where the conjectured holographic map could be defined. Trying to delineate this region is one 
of the things we attempt in the second part of this paper, although more work remains to be 
done. 

A very interesting feature of the numerical solution found by Choptuik [7] is that when 
p ^ p* the metric near r = is approximately discretely self-similar (DSS). This means that 
the metric functions satisfy the property Z^,{t,r) ^ Z^:{e^t,e^r), with A ^ 3.44 in the case 
of the spherical collapse of a massless scalar field in four dimensions. This is known as the 
"echoing" phenomenon. A similar analysis for the case of a perfect fluid collapse [S] shows that 
the critical solution exhibits continuous self-similarity (CSS) near the center. This implies that 
there exists a vector field v such that 

^vgf,u = ^fiVy + VyVf, = 2g^^, (1.2) 

where Cy represents the Lie derivative with respect to the vector field v. 

Space-times with CSS are very interesting from various points of view (see [H] for a com- 
prehensive rewiew). Although DSS is a remarkable phenomenon in gravity, it seems to be a 
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p>p* 




Black hole 

Figure 1: Phase space picture of the critical gravitational collapse. 

disadvantage when trying to establish a holographic correspondence with the Regge region in 
QCD, where no echo behavior is to be found. However, the so-called leading (log s) -behavior 
of the amplitudes does indeed show scale invariance (see Section H]). This is a fundamental 
reason to abandon the construction of the holographic map using collapsing massless scalar 
fields and to use instead a system where the critical solution exhibits scale invariance. The 
archetypical system of this kind is the spherical collapse of a perfect fluid. 

One of the main (technical) difficulties in the original computation of the Choptuik exponent 
[7] is that it requires a very involved numerical solution of the Einstein equations. In [HI [10] an 
alternative procedure to compute 7 was proposed based on a renormalization group analysis 
of critical gravitational collapse. In this picture, the surface p = p* represents a critical surface 
in the space of solutions separating the basins of attraction of two fixed points, corresponding 
respectively to Minkowski and the black hole space-times (see Fig. [1]). The critical solution 
with DSS or CSS has a single unstable direction normal to the critical surface. 

In this approach, the critical solution is characterized by having a single growing mode for 
perturbations around it. We can characterize it by the corresponding Lyapunov exponent. If 
Ai > represents the Lyapunov exponent associated to the repulsive direction, the critical 
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exponent 7 is given in terms of it by [TOl [6] 

7 = f. (1.3) 

Consequently, in order to compute the Choptuik exponent associated with a collapsing system, 
one should identify the critical solutions with either CSS or DSS and look for their single 
unstable mode. 

For systems in which the critical solution has DSS, this renormalization-group approach 
provides a nice geometrical picture to critical collapse. However, from a practical point of 
view, the computation of the DSS critical solution still remains a fairly formidable numerical 
task. If, on the other hand, we deal with systems exhibiting CSS the problem of finding the 
Choptuik exponent is reduced to numerically solving a system of ordinary differential equations 
and studying their linear perturbations. For this it suffices to use standard numerical routines 
to solve the corresponding system of ordinary differential equations, for example Runge-Kutta 
methods. This is a substantial simplification with respect to simulating the whole system using 
dynamical triangulations. This is the approach we follow in this paper. 

This paper is organized as follows: in Sections [2] and [3] we present the gravitational part of 
the problem. We start with the study of self-similar, spherically symmetric, solutions to the 
Einstein equations in d dimensions, sourced by a perfect fiuid with equation of state p = kp. 
These can be interpreted as critical solutions for the spherical gravitational collapse of a perfect 
fiuid. We write the corresponding equations of motion and study in detail the regularity 
conditions to be imposed on the solutions. Next, we integrate these equations numerically and 
study the properties of the solutions. This part has a value on its own. To our knowledge 
this is the first time such an analysis is carried out. We also present a number of technical 
clarifications that should be useful to those who would like to engage in similar computations. 

Before proceeding to numerically integrate the differential equations there are a number 
of non-trivial issues that need to be resolved. The Einstein equations with CSS have three 
singular points in terms of the scaling variable z. The solutions near two of them, z = 
and z = 00, can be dealt with using standard Frobenius analysis. The third, the so-called 
sonic point, is far more subtle. The generic solution present sonic booms where part of the 
collapsing matter is thrown out to infinity and part forms a compact object near the origin. 
For arbitrary initial conditions it is not possible to analytically join the solution below and 
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above the sonic surface. The critical solution is precisely the one that can cross analytically the 
sonic point (surface). This analyticity constraint implies that there is only a discrete number 
of possible solutions. The absence of "moduli" in the critical solution adds to the difficulty in 
its determination. The gravitational part of this paper is essentially the story of how to cross 
the sonic point, both for the critical solutions as well as their perturbations. This turns out to 
be a long and winding road we have tried to make as practicable as possible for the reader. 

Once the self-similar solution is obtained, we study in Section [3] its perturbations and 
compute the corresponding repulsive Lyapunov exponent for diverse values of k and d. From 
them we obtain the Choptuik exponents in each case, in particular the ones quoted above. 

After the detailed discussion of the critical gravitational collapse of a perfect fluid in higher 
dimensions, in Section H] we discuss its holographic interpretation in terms of four dimensional 
high energy scattering in QCD in the perturbative Regge limit. 

In Subsection 14 . 1 1 we provide a general introduction to the concept of the hard Pomeron as 
a bound state of two Reggeized gluons. The Reggeization of the gluon at high energies calls for 
an interpretation in terms of an effective string theory. We are always working in the small 't 
Hooft coupling limit in the gauge theory side. However, in deep inelastic scattering (DIS), the 
gluon distribution function is sensitive to nonperturbative effects not related to confinement: 
the so-called saturation effects. They are responsible for the transition from a dilute regime, 
where parton-parton correlations are weak, towards a dense limit where these correlations are 
very strong. We argue that this dilute-dense transition corresponds to the formation of a higher 
dimensional tiny black in the gravity side of our holographic proposal. 

In Subsection 14.21 a brief introduction to the generalities of evolution equations in DIS is 
provided. We then focus on the small-x limit and describe in some detail the building blocks 
of the BFKL (Balitsky-Fadin-Kuraev-Lipatov) formalism [14J. The BFKL equation describes 
a linear evolution of the parton densities with respect to x with no ordering in the virtualities 
of the internal propagators. This generates a UV/IR symmetric random walk in transverse 
momentum space which can be described by a simple diffusion equation. We identify the 
variable 1/x, proportional to the center of mass energy, as the holographic extra dimension. 

The saturated region, corresponding to a dense system of packed partons, is dominated by 
nonlinearities required to respect the Froissard-Martin bound on total cross sections for very 
large energies. In Subsection 14. 31 the simplest generalization of the BFKL equation containing a 
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nonlinear term, the BK (Balitsky-Kovchegov) equation [15j, is introduced. However, we do not 
use the BK equation in this work since it is posible to study the transition to saturation using 
an absorptive barrier in the BFKL equation in the hmit of asymptotically large center of mass 
energies. This barrier introduces an asymmetry in transverse momentum space suppressing 
the propagation towards IR virtualities. These correspond to gluons with a large transverse 
size and a larger probability of destructive interference. In Subsection 14.41 we discuss the 
holographic relation between saturation and critical gravitational collapse and argue that the 
suppression of IR modes in BFKL evolution corresponds to the large fraction of initial matter 
which evolves into flat space-time in the critical formation of the tiny black hole. 

To close the paper, in Section [5] we do not present conclusions, but rather a collection 
of unanswered questions and work in progress. For the reader's convenience some technical 
details have been deferred to the Appendices. 

2 Self-similar perfect fluid solutions in d dimensions 

This is a very long Section where we try to collect in a self-contained way our results on the 
gravitation collapse of perfect fluids in diverse dimensions. Although we are interested in the 
case of a conformal fluid in five dimensions, the complexity of the problem is similar in any 
dimension. Hence we have analyzed the problem including the space-time dimension c? as a 
parameter. The basic results are plotted in Figs. [2] to O where we exhibit the properties of 
the critical solutions for a sample of values of k in five dimensions. Similar plots could have 
been shown in other dimensions. 

In order to find the Choptuik exponent for the critical collapse of a perfect fluid in d 
dimensions we need to study the corresponding critical solution. In our case this is given by a 
spherically symmetric self-similar solution of the Einstein equations sourced by a perfect fluid. 
We begin therefore by obtaining the corresponding Einstein equations. 

It is important to stress that one of the difficulties in determining the critical solution in 
gravitational collapse is that it does not depend on any parameter; or in modern parlance, 
it generically has no moduli apart from the speed of sound and the dimension of space-time. 
As a consequence, we have to deal with a system of non-linear differential equations with 
three singular points that do not seem amenable to analytic treatment in spite of (our) many 
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unsuccessful efforts. In order to obtain numerical solutions that are regular at the singular 
points we have to toil through a number of steps that we have tried to make as palatable as 
possible. Once the critical solutions are determined we proceed to study their perturbations 
and to compute their Lyapunov exponents in the next Section. 



2.1 The Einstein equations 

Our starting point is then the spherically symmetric d-dimensional line element 
ds'^ = -a{t, rfdt^ + a{t, rfdr^ + R{t, rfdn] 



'd-2? 



(2.1) 



where dQ'^_2 is the line element on the {d — 2)-dimensional sphere S'^ ^. We look for solutions 
to the Einstein equations sourced by a perfect fluid with energy-momentum tensor 



T^,u = (p + p)u^,Uy + pg^y 



and equation of state 



0<k<l. 



(2.2) 



(2.3) 



We have written the line element (12. ip using comoving coordinates where the components of 
the fluid velocity are = (u*, ). This is quite standard in the literature. The basic advantage 
is that the integration of the Bianchi identities is straightforward. 



The equations of motion. The details of the calculation of the Einstein equations for the 
metric (12. ip can be found in Appendix A. As explained there it is enough to consider the 
(t, t), {t, r), (r, r) components given by 
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pR,tR''-\ 



(2.4) 



where Gn represents Newton's constant. The first equations implements the comoving condi- 
tion with vanishing space-like velocity components for the fluid. These equations have to be 
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supplemented by the Bianchi identities V^T^*^ 



a P + P 

^ = (2.5) 

a p+p ^ ' R ^ ' 

One of the basic simphfying features of studying spherically symmetric gravitational col- 
lapse is the fact that gravitational waves are not generated, hence we can write the metric 
components in terms of the energy-momentum explicitly. General collapse, without spherical 
symmetry is a far more difficult problem that has not been analyzed yet in the search of critical 
phenomena. As in the four- dimensional case [llj, the equations in (12.41) can be rewritten in a 
simpler form by introducing the function m(t, r) defined by 

2GNm = R''-''(l + ^-?^\. (2.6) 

Together with the Bianchi identities (12. 5p we arrive at the following system of equations 

X StT d, 2i 

'^^,r = T; Gn pR,rR , 

a — 2 

d-2 



2m,i 


IGvr 

= ( 

d-2 


2Gis[m 






P,r 


a 


p + p' 


a,t 


P,t 


a 


P + P 



Gjy pR fR 

1 + ^-%), (2.7) 



These equations plus the equations of state (12. 3p make a total of six equations to determine 
the five unknown functions, a{t,r), a{t,r), R{t,r), p{t,r) and p{t,r). There are therefore only 
five independent equations. Their physical meaning is transparent. The first one indicates how 
the mass function m defined in Eq. (12.61) grows with the positive energy density p. The second 
implies the decrease of this same mass function induced by a positive pressure. The third 
gives a measure of the formation of an apparent horizon. Notice that equation (12.61) can be 
written as the condition when the surface R{t, r) = constant become null. In these coordinates 
it represents an integral of motion. 
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The Bianchi identities can actually be written as 

k 



(loga),t = -■^(logp),t-(rf-2)(logi?),„ (2.8) 



hence 



ait, r) = c„(t)p-^, a(t, r) = Ca{r)p-^^ R^~'', (2.9) 

where Ca{t), Ca{r) are arbitrary functions of t and r respectively. The function Ca{r) can be 
absorbed by a redefinition of the radial coordinate r. Thus we set Ca(r) = 1 in the following. 

Special velocities. There are two families of surfaces in the geometry given by the line 
element ( 12.11) . The first ones are 

i?(t,r) = constant, (2.10) 

defining spheres of constant area. The second family of surfaces comes about because we are 
interested eventually in self-similar solutions. They are defined by 

r = —zt, z = constant. (2-11) 

One can define now the velocity of these surfaces as the velocity of a fiducial observer lying on 
one of them with respect to the fiuid. This is given by the ratio between the proper distance 
di = adr that this observer covers and the proper time d\ = adt it takes. In the case of the 
constant-i? surface we find 



= ^dX + ^di ^ Vn = -^, (2.12) 
a a aR, 



whereas for surfaces of constant z the velocity is 



= --dX + -di, =^ V, = - — . (2.13) 
a a a 
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Changing coordinates. Since eventually we are interested in self-similar solutions it is 
convenient to define the new coordinates (r, z) by 

r=-log(-l). (2.14) 

where we have introduced a length scale ig in order to make the argument of the logarithm 
dimensionless. Notice that the collapse takes place for t < 0, the critical black hole being 
formed in the limit t 0^. This corresponds to the limit r oo. In addition, we introduce 
the dimensionless functions rj{t,r), S{t,r) and M(t,r) defined respectively by [11] 

SiT,z) = (2.15) 
r 



M{t,r) 



2GNm{t,r) 



Up to this point we have written the Newton constant Gn explicitly. From now on we use 
units in which G^v = 1- In this system of units, p and p have dimensions of (length)^, whereas 
m is measured in units of (length)*^"^. The powers of the Newton constant can be restored in 
all expressions by replacing m G]\fm, p G^p and p Gj\fp. 

Denoting now derivatives with respect to r and log z respectively by a dot and a prime, 

d , d 



dr ' d log z ' 

the Einstein equations become 



(2.16) 



Sd-3 




M 



where the dimensionless function ?/(r, z) is defined by 



y = -{d-2){d-3)^^ (2.18) 
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and Vz is given by Eq. fl2.13p . Both a(t,r) and a(t,r) have been integrated in fl2.9p and can 
now be expressed in terms of the new functions as 

a = CaiT)ljy^\ a = r]-T^^S^-''. (2.19) 

Thus the velocity is known and we are left with three functions M(r, z), S{t, z) and y(r, z) 
[or ?7(r, z)] to be determined from the three differential equations in (I2.17p . 

Self-similar solutions. We now impose self-similarity in the metric (12.11) with respect to 
the conformal Killing vector field 

v = tdt + rdr (2.20) 

that satisfies Eq. (11.21) . In the new coordinates (I2.14p the action of the homothety generated 
by V simply shifts the new time variable r r + constant. Thus, the condition of self-similarity 
translates into the condition that all dimensionless functions in the problem are independent 
functions of r and depend only of the scaling coordinate z. This is the case for the functions 
defined in Eqs. (I2.15p and (I2.18p . Then, the Einstein equations are (cf. pT] ) 

_ id-3)k S' 
M ~ y S' ^ ^ 

Using the first two equations we can solve for M' and S'. The third one can be interpreted as 
an equation defining ri{z) as a function of S{z) and M{z). Taking its derivative with respect to 
log z and after a tedious computation one arrives at the following system of Einstein equations 
dlogM {d-3)k (I ^ 



d\ogz k + 1 \y 

dlogS _ 1 
d log z k + 1 

d log 7] 1 

d log z — 



(y-l), (2.22) 
[l + kf ^ 



d-2 

12 



T]^S^-'^ -{d-2){y- 1)V; - 2k 



Since in deriving the last equation we have taken a derivative these equations have to supple- 
mented by the third equation in (I2.2ip at least at a particular point, and also by the expressions 
for a{z), a{z) and Vz{z) 

( z^ \ ^ 

«W = ' (2-23) 

a{z) = r^-^S'-^ 

V.iz) = --.'^ (^) ^'-^ 

where we have fixed the gauge by setting Co(r) to a constant Ca- Notice that the system 
still has a residual gauge freedom corresponding to constant rescalings of the time coordinate 
t (t, with C 7^ 0. This implies a rescaling of the constant by — > (Ca- Without loss of 
generality we require a{0) = 1. A very important equation is the form of the velocity field Vr 
for self-similar solutions. From fl2.12p and the self-similar ansatz we can write 

Vn = -V. = -V; 1^. (2.24) 

This equation implies in particular that, if Vz 7^ 0, the zeroes of Vr are in one-to-one corre- 
spondence with the values of z where y = 1. 

The sonic point. The Einstein equations (12.221) present a number of singularities. In princi- 
ple the possible singularities at z = and z = 00 are expected from the fact that we are looking 
for self-similar solutions. As in the case of regular singularities of linear differential equations, 
these singularities imply that the solutions have a Frobenius behavior in their neighborhood. 
Hence they behave as z"'f{z), where the power a depends on whether we are at 2; = or 2; = 00, 
but f{z) is in any case an analytic function. 

The singularity at = k, the sonic point or sonic surface, is more delicate. It represents 
the locus where the velocity of the lines of constant z reach the speed of sound. At this place 
shock waves can form. This means that if we start integrating the solution from the origin (or 
from infinity), we may not be able to cross the sonic surface analytically. Geometrically, the 
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line z = Zsn where Vz( Cauchy horizon. Choptuik's critical solution is everywhere 

analytic, hence we need to impose stringent conditions at the sonic point so that the solution 
passes through smoothly. These conditions will be explained in detail later. 

A generic solution of the system (12.221) depends on three initial conditions. When requiring 
analyticity throughout, we start by imposing it at the origin. As we will show below this 
reduces the family of solutions to a one-parameter family. If we also require analyticity at 
the sonic point, we find that only a discrete set of values of this parameter yield admissible 
solutions. Generically, the critical solutions have no moduli. They are classified by the zeroes 
of the velocity field (12.121) . In particular, the critical solution we are interested in, the so- 
called Hunter-A solution |13J, is characterized by having a single zero of Vr p!2| [TT] . The lack 
of moduli characterizing critical behavior is one of the main difficulties in trying an analytic 
approach to the solutions of interest. We can easily compute their asymptotic behavior, but 
close to the sonic point, unfortunately, they only seem to be amenable to numerical analysis. 

It is reasonably easy to compute the asymptotic properties of the solution at the origin and 
at infinity. However if we want to study the solution near the sonic point, we need to impose 
some restrictions there. In particular regularity of the third equation in (12.221) at z = Zgp 
imposes the additional constraint 



d-2 



0. 



(2.25) 



This, however, is not enough to determine (logr^)' at the sonic point, since the right-hand side 
of the corresponding equation is still indeterminate. In order to evaluate this derivative we 
have to define the right-hand side at the sonic point using I'Hopital rule 



dlogT] 



dlogz 



1 



d-2 



77^5^-2^ -(rf-2)(y-l)lf 



2k 



(2.26) 



As we will show in detail in Subsection 12.21 evaluating the derivatives on the right-hand 
side of Eq. (12.261) and using the other equations of motion leads to a quadratic equation to 
be satisfied by (log?])' at the sonic point. This means that one has to determine which of the 
two solutions provides the correct analytic solutions connecting smoothly the z = to z = oc 
through the sonic point. The "wrong" solution would show catastrophic implosion to form a 
black hole, an explosion to shed all collapsing matter to infinity or something in-between. 
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In the critical solutions, on the other hand, there is a balance between the matter that 
collapses and the matter that is ejected, so that eventually we only form a black hole of 
"zero" mass. It can be shown that once the first derivative of (log?7)gp is determined, all 
higher derivatives can be calculated unambiguously at the sonic point through linear equations 
obtained by repeated application of the I'Hopital rule. Indeed, the function logr] can be 
expanded near z = Zgp in a power series in \og{z/zsp) 



where i/sp = y{zsp). The coefficients (^/gp) are calculable from fl2.26p and the equations 
obtained by taking further derivatives with respect to logz in the third equation in (12.221) . 
In Subsection 12.21 we evaluate the first coefficient of the expansion, Ci^^\ysp) = (log?7)sp. The 
details of the calculation of Cij^\ysp) = (log?7)gp are given in Appendix B. 

From the previous discussion it follows that one of the main problems we have to solve is 
to determine the correct sign in the quadratic equation for (log?7)sp. This input is crucial in 
the numerical evaluation of the critical solution. Before we get there, however, we need to get 
acquainted with some asymptotic properties of the solutions to Eqs. (I2.22p . 

Behavior near z = and z = oc. A good place to start the numerical integration of (I2.22p 
in the search of the critical solution is the origin (or infinity). Hence it is worth studying our 
system of differential equations in these regions. Alternatively, we could start the numerical 
integration at the sonic point and integrate towards z = and z = oo. Although there are 
some advantages in doing this, we stick here to z = and defer the analysis of this other case 
to a future work [T6] . 

We start the analysis by finding the conditions imposed by analyticity at z = 0. By 
demanding regularity at the origin, we find that the three-parameter family of general solutions 
to (I2.22P reduce to a one-parameter family. 

The asymptotic behavior as z ^ can be found from the first equation in (12. 7p . Since 
for a regular solution R{tQ, r) is expected to be a single valued function of r, Eq. (12. 4p can be 
formally integrated at fixed time t = to < 




(2.27) 




(2.28) 
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For small values of the r-coordinate, r 0"'', one expects the coordinate R{t,r) to vanish. 
Then the integral can be approximated by 

87rp(to,0) 



m{to, r ~ 0) 



-i?(to,r ~0) 



d-l 



2.29) 



{d-l){d-2) 

Writing now the function y{t,r) defined in Eq. (2.18) in terms of m{t,r), p{t,r) and R{t,r 
we find its limit when 2; ^ 0^ to be 

y{t,r) = {d-2){d-3) 



m 



d-3 



(2.30) 



SirpRd-^ ' d-l 

It is important to notice that this behavior of the function y{t, r) as r — > is completely general. 
It does not require the solution to be self-similar, and it is valid as well for the perturbations 
to the self-similar solution. It only depends on the regularity of the solution at the origin. 

Once the limiting value of y{z) as z ^ 0^ has been obtained it is possible to solve the 
first two equations in the system (12.221) for small z. The fact that the value y{0'^) is universal, 
together with the quadratic relation (]2.2ip . leads to the fact that the asymptotic behavior of 
the functions M{z), S{z) and ri{z) depends only on a single free parameter. Following [11], we 
use the constant D defined by 



2D = lim 



8nt^p{t,r) 



(2.31) 



From Eqs. (12.221) in the limit z ^ 0^, together with the gauge condition on the metric 
coefficient a{z) — ^ 1 in the same limit, one finds the asymptotic form of M{z), S{z) and t]{z) 



near z ^ 0^ as 
M{z) 1 

S{z) : 

ri{z) : 



^ I (2D)T^ 

Zk+i 



id -2) 



k+1 



{d-l)k + d-3 



Z (d-i)(fc+i) 



{2D)T^ 



k + 1 



/, d-3 

k + 



d- 1 



n=l 



l-d 00 



(n) o„ d-3+fc(d-l) 



(2.32) 



n=l 



2D + j2vl; 



^ '2 (d-l){fc+l) 



n=l 



Apart from an overall power of z, all functions admit an expansion in integer powers of the 

^ d-3+fc(d-l) 

variable z (d-i){fe+i) . It should be noticed that the exponent of the expansion variable never 
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vanishes for physical values of k. Incidentally, the constant Cq, can be written in terms of 

fc 

D as Ca = (2-D) . It is clear from this expression that D is a gauge dependent quantity 
with respect to the residual gauge transformations z X^^z, Ac^. We conclude from 

this exercise that analyticity at the origin gives a family of solutions depending on a single 
parameter D. 

A further restriction on the possible values of D comes from requiring that the solution is 
analytic at the other singular value of z, the sonic point z = Zgp. Generically we find that only 
a discrete set of values of D lead to a regular solution for all finite values of z. These correspond 
to the higher- dimensional generalization of fiat-Friedmann, general relativistic Larson-Penston 
and Hunter-A/B solutions [11]. How these solutions are found numerically is discussed in detail 
in Sec. (ESD- 

The behavior of the functions in the limit z ^ oo can also be obtained using similar 
techniques. It can be shown that all functions M{z), S{z) and ri[z) tend to constants in that 
limit and that from the equations of motion f l2.22p it follows that y{z) — 1 when z oo. 
Looking at the subleading terms, the expansion is of the form 

Z{z) = Z^ + Z^^^z-"^^+..., (2.33) 

for any of the relevant functions collectively represented here by Z{z). All functions near 

i-fc 

z = oo can be expanded in inverse power series in the variable z . This expansion parameter 
is different from the one appearing when expanding the functions around z = 0. It is important 
to notice that this expansion parameter is independent of the dimension. Furthermore, if we 
look at the expression for in (2.23), we find that its absolute value diverges at 2; = oo 
precisely as ^1+*= . In fact, V/ provides a good uniformizing variable for the system of equations 
under study. Indeed, is a monotonic function in the whole region between z = and z = 00. 
This variable is essentially the one used in the analysis of Ref . . We turn now to study the 
form of the field equations in terms of this variable. This will be helpful in determining the 
correct sign of (log?])' at the sonic point. 

The uniformization variable. One unpleasant aspect of the equations obtained so far is the 
presence of generically irrational exponents in both the equations and the initial conditions. 
This problem can be solved by choosing an appropriate uniformization variable in terms of 
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which integer powers are obtained. This can be done by noticing that the expansion of V"^ 
around z = has the form 

2^ 



V 



r, d-3+k(d-l) 
2; (d-l)(fc + l) 



(2D)-7 



d-2)(fc+l) 



d- 2 



) d-2 
'd-1 



(2.34) 



n=l 



We have seen before how all functions, apart from an overall power of 2;, can be written as 

2 d-3+fc(ii-l) 

series expansion in z (d-i){fc+i) _ Therefore, by redefining the functions properly in such a way 
that all the Frobenius factors are absorbed, it should be possible to write them as a series 
expansion in integer powers of V'^ . Moreover, for critical solutions Vzlz)"^ is a monotonous 
function of z that can be used as the independent variable. 
We begin by defining a new function S{z) as 

S = -{2D)-T^z'^V,S^~^'^. (2.35) 

Using this new function together with y{z) and Vz{z), the system ( 12.22^ can be rewritten as 



dhgS 
dlogz 



dlogy 
dlogz 



dlogz 



2{d-2) 
k + 1 



k{d-3) (I 



2k+{d- 2)V^{y - 1) - 2D ^^^^}^ S 



d-2 



k+l 



+ 



y 



jd-l) 
k + l 



2 -2k 2{d-2 
1 



2k + {d ~ 2)V^{y - I) - 2D 



{k + lf 
d-2 



S 



(2.36) 



iy 



f 1 k + l 

2{k-l) 
{k + l){V^-k) L 



2k + {d - 2)V^{y - I) - 2D 



{k+l)- 
d-2 



-S 



After the introduction of S{z) all functions appear with integer exponents on the right-hand 
side of the equations. At the same time the parameter D appears explicitly in the system to 
be solved. 

The velocity Vz{zy' vanishes in the limit 2; — 0^. Similarly, the function S has a simple 
behavior in that limit function of V^, namely 



(2.37) 
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This shows that is a good candidate to be used as independent variable. With this choice, 
the functions ^(V^^) and y{V^) satisfy the following system of differential equations 



dhgS 



-2(rf-2)(y-l)-^ 


'2k + {d-2)V,'{y-l)-2D^S 




2~2k-2id-2)iy-l)-'-^ 


'2k + id-2)V,^{y-l)-2D^S 



(2.38) 



dlogy 


k{d-3) (i-l) 


-{d- 




1 \ I fc+i 
^) + V,^-k 


'2k + {d-2)V,\y-l)-2D^S 


diogv;2 


2 - 2A; - 2{d 


-2)(y 


-1)- 


2{k-l) 
Vi~k 


'2k + {d-2)V^{y-l)-2D^^S 





This system can be integrated in principle by requiring regularity at the sonic point = k. 
Notice that in these variables the regularity at these points leads again to the condition (12.251) . 
As a bonus, using these variables we eliminate the constant D from the initial conditions at 

= to make it appear explicitly in the equations. In integrating the system numerically we 
can use either Eqs. (I2.22p or Eqs. (12.381) . Both systems are equivalent, although the conceptual 
advantage of the parametrization is that near the regular singular points the functions can 
be expanded in integer powers of this variable. 

Once the functions S{V^) and y{V^) are known, it is possible to find their dependence on 
the scaling coordinate z by integrating the last equation in (12.361) to find Vz = Vz{z). This 
procedure gives the same values of D for the regular solutions as the ones obtained integrating 
the equations in terms of the independent variable z. The fact that the values of D are the 
same is a consequence of the fact that in writing the equations (I2.38P we have not changed the 
gauge fixing conditions. 



Shooting from the sonic pointo. The system of ordinary differential equations (12.221) can 
also be integrated from the sonic point z = Zgp. In this case it is convenient to fix the residual 
gauge freedom, t (t, by locating the position of the sonic point at a given value of the 
scaling coordinate, 2;sp = ^ > 0. In this case the constant Cq, is fixed by the equation defining 

k 

r, 1-fc 



the sonic point, V^(-2sp)^ 



i+k q4— 2d 

sp 5 



(2.39) 



^This paragraph can be read independently of the main argument hne. The reader can skip this paragraph 
without losing the thread of the arguments. In a future publication [16j we plan to present the numerical 
integration from the sonic point towards the origin and infinity. 
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where the subscript "sp" indicate that the corresponding function is evaluated at the sonic 
point. 

In the same way as the behavior of the solutions near the origin is determined by a single 
parameter D, the requirement of regularity at the sonic point forces now that the values of the 
functions at this point depend on just one free parameter. We can take this to be the value of 
the function y{z) at the sonic point, ?/sp- To see this we notice that there are two conditions 
that the functions should satisfy at the sonic point, namely the regularity condition (12.251) and 
the quadratic equation in fl2.2ip 

2k + {d-2)kiy,p-l)-^^j^r],Xp = 0, 

i^+krU-^,)+k-^ii'ys,r'-^ik+ys,r = o. (2.40) 

V '-'sp / ""sp '^sp 

The quantity Mgp can be now written in terms of r^gp by using the definition of the function 

1 

y{z) in Eq. fl2.18p . whereas = r/sp*^^ 5's In this way it is possible to solve the previous 
two equations for rj^p and S^p as a function of ysp. Moreover, by requiring analyticity of the 
functions at Zgp it is possible to show that all derivatives of the functions at that point are also 
determined by ysp. 

Thus, we have found then that there is a one-parameter family of solution to the system 
(12.221) regular at the sonic point. Only when we require in addition that the corresponding 
functions also show the behavior described previously as 2; ^ 0+ the set of allowed functions 
restrict to a finite set of values of the initial condition y^p. These allowed values for each 
k correspond to different types of solutions associated with the values of D giving regular 
solutions everywhere. 

It is interesting to relate the initial condition y^p at the sonic point with the parameter D 
governing the asymptotic behavior of the regular solution at 2; = 0. The key ingredient to do 
this is to notice that the condition Cq = (2D)'^ together with Eq. (12.391) gives the following 
relation between D and ysp 

k+l 

2^ = e^(^^srj (2-41) 

where as shown above Ssp and r^sp are given in terms of ysp. The explicit dependence of this 
expression on the "gauge fixing parameter" ^ shows once more that Z) is a gauge-dependent 
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quantity. In addition, from this expression it is immediate to verify that the combination 

fc-i 

D Zsp is gauge invariant. 

2.2 Dealing with the sign ambiguity. 

Let us go back to the determination of (log//)' at the sonic point. The last step that we need 
to solve before embarking in the numerical solution of the self-similar Einstein equations is to 
fix somehow the ambiguity associated with the choice of the right solution to the quadratic 
equation satisfied by (log?7)gp. This leads to a long and involved calculation. This is, in fact, 
one of the prices one has to pay for using comoving coordinates. In Schwarzschild coordinates 
a similar problem appears [10] but it is technically easier to solve. If the reader is willing to 
give us the credit that we have solved the problem, this subsection can be skipped. The result 
however is very important in the numerical determination of the critical solution. 

To see how the quadratic equation for (log?7)gp comes about we go back to Eq. (I2.26p . In 
order to evaluate the right-hand side of this equation we need the expressions for iV^y, {a^rf)' 
and y' evaluated at the sonic point. For the first one we use the expression of in Eq. fl2.23p 
to write 



2y 

z /sp 



l + k 1 + kri ^ ^ S 



V 



2k 



__[l_fc+(2-d)(y,p-i: 



mogvYsp] 



(2.42) 



where in the second line we have used the field equations fl2.2ip to write the logarithmic 
derivatives in terms of the values of the functions at the sonic point. For rja'^ we find also from 
Eq. (ESS]) 

1 



{r]a 



2\r 



/sp 



k + 1 



'7sp'^sp 



(fc-l)(logr^)' +2(2-rf)(ysp-l) 



(2.43) 



Finally, the derivative of the function y can be written as 

k{d-3) f 1 A d-1 



(logy)s 



sp 



1 



sp 



(log^)s 



sp 



(2.44) 



k+1 V?/sp / k + 1 
Substituting these expressions in Eq. (12.261) leads to the following quadratic equation for 
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sp 



2k r 
__ [1 _ + (2 - d)(y3p - 1) - (1 - fc)(logr/)y [(logr;)^p + {d - 2){y,^ - 1) 

k + 1 



+ 



id -2) 

{d-2)k 



k + l 



(Z/sp 



{k-l){\ogr,)'+2{2-d){y.. 



(2.45) 



(rf - 1) (i/sp + fc) - 2fcJ +{d- 2)A;?/,p(logr/)^p. 

As discussed above, the values of all functions at the sonic point can be expressed in terms 
of a single parameter that we can choose to be ysp. This means that the coefficients of the 
quadratic equation fl2.45p can be written in terms of this single parameter. To do this we solve 
Eqs. f l2.40p for the combinations a'^pTisp and a^pS^ with the result 

k{d-2) 



^spC^sp 



ik + l] 



{d-2)y,p-{d-4) 



[1 + kYaiS 



{k + Vsp? - - Vspf + 



2k 



"sp'-'sp V" ' fspy '"V^ fspy i ^ gZ/sp 

Substituting now these expressions in Eq. fl2.45p we arrive at 



{d-2)y,p-{d-A) . (2.46) 



k + {2-d){y, 



sp 



1 - k){\ogii)' (logr^)' +{d- 2){y,p - 1) 



[d - 2)y,p - (rf - 4) (A; - l)(log77)' + 2(2 - d){y,p - 1 



+ id-2) (ysp - 1) 



(2.47) 



{d - 1) (t/sp + k)- 2A;J +{d-2){k + l)|/sp(logr7)^p. 

As advertised, all coefficients depend only on ?/sp- 

This quadratic equation can be written in the canonical form as 

[(logr/)y2 + Aiysp){\ogvy.p + B{y,p) = 0, (2.48) 

where the coefficients are given by 

2{k-l)Aiysp) = M-2-{d + 2)k-A{d-2)y,p. 

2{k-l)B{y,p) = (rf-2)(ysp-l)[6-(d-l)(ysp + fc)]. (2.49) 

The first thing to be noticed is that the equation for (log?7)gp is independent of ^ and therefore 
it is gauge invariant, as it is to be expected. Solving the equation 

^(Z/sp) 



sp 



± 



4 



(2.50) 
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There is therefore an ambiguity in the value of (log?7)sp given by the two possible signs in this 
equation. 

The numerical integration of the equations (12.261) from the origin z = shows that, for both 
the Hunter-A/B solutions, (log?7)gp is determined by the minus branch in Eq. (12.501) . For these 
types of solutions, this branch is automatically selected by the system when regularity at both 
the origin and the sonic point is imposed, i.e. when the constant D takes the corresponding 
critical values. The situation changes if we start the integration from the sonic point. In this 
case we find that the "correct" branch in Eq. (12.501) has to be chosen in order to determine 
the right initial conditions for the numerical integration. This situation will be addressed in 
detail in a separate publication [T6] . 

2.3 Numerical analysis 

After this long study of the properties of the equations of motion (12.221) for the self-similar 
background, we can proceed to integrate the equations numerically. To do so we use a fourth- 
order Runge-Kutta routine that integrates the system starting at z = 0. As explained before, 
the main difficulty lies in finding the right value of D providing an analytic solution at the 
sonic point. For a generic value of D the routine gives a solution that diverges at this point. 

In order to find the right value we perform a scan in D to find the smooth numerical 
solution at the sonic point ^sp- Since we look for solutions that can be interpreted as describing 
critical gravitational collapse of a perfect fiuid, we are interested in the solution with a single 
unstable mode. In dimension larger than four, this corresponds to the higher-dimensional 
generalization of the general relativistic Hunter-A solution, which are characterized by single 
zero of the function Vr. 

The numerical analysis is a bit delicate due to the singular character of the sonic point, 
where minor numerical noise can easily lead one astray. Here it is important to use the 
analytical information we have at hand: the behavior of the solution at the origin and the 
value of the derivatives of the relevant functions at the sonic point. After long numerical work, 
we have determined the critical solutions for different values of k in various dimensions, in 
particular d = 4,5,6 and 7. We have included the d = 4 computation in our analysis in order 
to compare with previous work, in particular Ref. |10j . 

In Figs. [2] and [3] we have plotted the profile of the main functions for three values of k 
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in five dimensions. The position of tlie sonic point is denoted by a dot on tlie corresponding 
curve. We see tfiat tlie curves fiave a structure completely similar to their four-dimensional 
counterparts, as found for instance in Ref. [11]. 

Fig. m shows how the function V"^ is monotonic in the whole range of values of z. This 
demonstrates that this function can be used as a good independent variable to describe the 
evolution of the system. 

Finally, in Figs. [5] and [6] we plot the various functions of interest against y instead of logz. 
This is useful to show that we have indeed found the critical solution. As explained before the 
critical solutions are characterized by the zeroes of the velocity field Vr (see Fig. 3). What 
Figs. [5] and [6] show is a form of critical damping. The solution tends to its asymptotic value 
qX z = oo after turning around it only once. Critical solutions with more than one zero in 
would spiral around their value at infinity a number of times given by the number of zeroes of 
this velocity field. Notice that y cannot be used to uniformize the problem, because it is not 
monotonous as a function of z. However, the fact that y[^z) only crosses the value y = 1 (at 
which = 0) at a single finite value of z before monotonously decreasing to ?/ = 1 as 2; oo 
is also a sign that the solution is critical. Indeed, the relation between Vr and y given in Eq. 
(12.241) shows that for the critical solution can be also characterized as the one for which the 
function y(^z^ crosses the value y = 1 at a single finite value of z. 
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Figure 2: Plot of M{z) and S{z) for k = 0.01, k = 0.1 and k = 0.2 in five dimensions. The 
dot on the curves indicates the position of the sonic point. 
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Figure 3: Plot of r]{z) and y{z) for k — 0.01, k — 0.1 and k — 0.2 in five dimensions. The dot 
on the curves indicates the position of the sonic point. 
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Figure 5: Plot of M and S versus y for k = 0.01, A; = 0.1 and k = 0.2 in five dimensions. 
The circular dot indicates on each curve the position of the sonic point. The corresponding 
coordinates indicate the position at which the system is attracted as z — > oo. 
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Figure 6: Plot of ri and versus y for k = 0.01, A; = 0.1 and k = 0.2 in five dimensions. 
The circular dot indicates on each curve the position of the sonic point. The corresponding 
coordinates indicate the position at which the system is attracted as z — > oo. 
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3 Perturbations of the self-similar solutions 



Once we have obtained the critical self-similar solutions we proceed to study their stability by 
considering perturbations around them. The perturbations break self-similarity and therefore 
have an explicit r dependence. In the context of critical gravitational collapse we are interested 
in solutions with a single growing mode. The analysis is long, and in particular we need to use a 
very relevant equation that does not seem to appear in the literature (see the paragraph on the 
?7-equation) . The sonic point, together with a new singularity that appears in the perturbation 
equations at y{z) = 1, play an important role in the analysis. Similarly to the calculation of 
the self-similar brackground, analyticity at all singular points gives the key to the computation 
of critical exponents. 

3.1 Analysis of the perturbations. 

The equations for the perturbations. Let us denote by Z{t, z) any of the functions 
involved in the problem. We consider perturbations around the self-similar solution Z{z) of 
the form 

Z(r, z) = Z{z) [1 + e e^" + ...], (3.1) 

where A is the corresponding Lyapunov exponent. Since we are interested in the mode that 
grows as r — > oo, the unstable mode is associated with a positive value of A. The dots in the 
equation include the contribution of stable modes with negative Lyapunov exponents. 

In order to write a set of differential equations for the perturbations Mi{z), Si{z) and rii{z) 
we have to go back to the full equations of motion 02.221) and plug in the ansatz fl3.ll) . Then 
we collect all the terms linear in the expansion parameter e. From the first two equations in 
(12.171) we find the equations of the perturbations Mi{z) and Si{z) to be 

{l + k)yM[ = {3-d)kyi~X[Miy + {d-3)kSi], 
il + k)S[ = yy,-^[M,y + id-3)kS^]. (3.2) 

We still lack one more equation for the perturbation of a third function that we chose to be 
yi{z). To find it we have to use the third equation in (12.17^ and expand to first order in e. 
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After a long calculation one arrives at 



Qyi = VmM^ + VsSi, (3.3) 
where all the "coefficients" are functions of z given by 

Q = {V^-k){l-y){k + y), 

+ ^^^^y[V^l-y) + k + y], (3.4) 

Vs = [l-k{d-2)]V,\l~yr-kid-l){k + yr+^-^^yr^i^^S'-'' 

+ X{k + 1) [V^\y-l) + k{k + y)]. 

Expressions ( 13.21) and (13. 3p give the relevant equations governing the perturbations. We 
can actually reduce them to a system of two differential equations by plugging the value of the 
perturbation yi, obtained from Eq. (13. 3p . into Eqs. (13.21) . As a result one gets the following 
system of two first order differential equations for the perturbations Mi{z) and Si{z) 

S[ = AsMi + BsSi, (3.5) 

where Am, 13 m, As and Bs are given respectively by 

{d-3)kVM + XQy 



A 



M 



Bm = -{d-3)k 



{k + l)Qy 
Vs + \Q 



{k + l)Qy' 



{d-3)VM-XQ 
= \d-3){k + l)Q^ ^^-'^ 



B 



s 



yVs -XkQ 
{k + l)Q • 
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Singularities and the »7-equation. As in the case of the background equations one has to 
be careful about potential singularities in the equations for the perturbations. In particular, 
by looking at the expression of yi in Eq. (13. 3p and the form of the coefficients Q, T>m and P5, 
we find that there are singular points corresponding to the zeroes of Q. This can happen at 
three points. The first one is again the sonic point z^p where = k. The second point is the 
value of Za at which y{za) = 1. The third possibility, y{z) = —k, is never realized for physical 
solutions with k > because the function y{z) is always positive definite. 

In the case of the self-similar background, regularity at the sonic point implied the extra 
condition fl2.25p to be satisfied by the physically admissible solutions. In the case of the 
perturbations an extra condition can be also extracted by considering the perturbations of the 
quadratic equation in (12.211) in a consistent way. In particular we take the derivative of this 
equation with respect to r and express the result in terms of the functions S{t, z), r]{T, z) and 
y(r, z). After these manipulations the resulting expression can be written as a linear equation 
for T] of the form 



where the coefficients Ci and C2 are given in terms of the functions M{t,z), S{t,z) and their 
derivatives. 

For the self-similar ansatz, Eq. (13.70 reduces to the second equation in (12.221) . If, on 
the other hand, we consider perturbations of the self-similar solution this expression contains 
important information about the regularity of the perturbed solution at the points at which Q 
vanishes. To see this, we insert the perturbations (13.11) into (13. 7p and use Eqs. (13.20 and (13. 3p 
to obtain a relation between Mi and Si. A long computation yields 



CiT] + C2 = 0, 



(3.7) 



GmMi + GsSi 







(3.8) 



where Qm and Qs are respectively given by 



Gm 



2{k + l) 
d-2 



k+ll 



2A r 



V,' + k-id-3)id-2){V,'-l)y As 



(3.9) 



+ 



(k + l) 



2A 




^ {(1 - k)V^ [{d - 3)(d - 2)y -lf-{k + ifS'-^'v-^^ {2S'vy - l) } 
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and 
Gs 



+ 



+ 



2A 

k + 1 

2A 

k + 1 

2A 



{k + l)Bs + {d-3){d-2)y-l 



S' 



4-2cZ 



+ 



2A 

k + 1 



Vs-{d-l)kQ 



Q 



fjk+i 



k{k + l) + {d- 2){3kd + d - 5k - 3)y - 2{d - 2)-^y 



s 



i-2d 



(3.10) 



V^ + k-{d-3){d-2){V!-l)y 



2A 

k + 1 



V 



{d-2){d-3)y 



2dk-3k-l + {k-l)--f 



id-3){d-2)y 



(k + l)^-' 

The coefficients Q, T>m, T^s^ and As are the ones defined above. We see that the only 
poles in Qs and Qm come from the zeros of Q, i.e. the points at which y = y^p, y = 1 and 
y = —k. Thus, Eq. (13. 8p can be written as a sum of terms potentially divergent at Q = plus 
those which are regular at those points. After some manipulations we find that this equation 
has a structure of the form 



^5iiML±^ + ,egula.. terms = 0. 

where we have introduced the coefficient defined by 



d-2 
2A 



yrjk+i b 



A-2d 



2A 

{k + iy^ 



V^ + k-{d-3){d-2){V^-l)y 



(3.11) 



(3.12) 



— |(A; - 1)V^ [{d -2){d- 3)y - l] ' + (A; + ifS'-'^r^-^^ {2S%v - l) } • 



{k + 

In order to preserve the identity (13. lip one must require that the whole expression is regular 
at the zeroes of Q. Since numerically it can be shown that JF never vanishes at those points 
we arrive at the conclusion that 



0. 



(3.13) 



zeroes of Q 

In particular, this equation has to be satisfied at the sonic point and provides a nontrivial 
condition to be satisfied by the perturbations there. 

Notice that this condition has been derived here as a straightforward consequence of the 
perturbations of Eq. (13.70 . which becomes trivial after imposing the self-similar ansatz. It 
should be said that Eq. (I3.13P plays a basic role in the determination of the Lyapunov exponents 
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of the critical solutions, and that to our knowledge it has never appeared in previous literature. 
Without including it, the path to understanding the properties of critical perturbations is rather 
treacherous. A systematic computation of the Choptuik exponents for different values of k and 
d relies heavily on it. 



Behavior of the perturbations at the sonic point. In order to negociate the singularities 
of the differential equation (13. 3p we impose analyticity at the singular points. For this we have 
to use again the I'Hopital rule to define the perturbation yi{z) in terms of the corresponding 
ones for the other functions. This means that in the regions near the points where Q{z) = 
we can replace Eq. (13. 3p by 



Q'y, = {VMM, + VsS,y. 
At the same time, instead of using Eq. (13.51) we consider the modified system 

S[ = AsM, + BsS^. 
The new coefficients Am, Bm, As and Bs are given by 



(d-3)fc , 



y 

{d-3)k 

y 



V's + \ 



TZBm = 

TZAs = y^V',i-\ 

UBs = yV's-Xk 
where the function TZ{z) is defined by 



k + y 
k + 1 



Q> k + y ■ 



d-3 ik + l)y 



{k + l)y 



k 



7^ = (A; + 1)Q' - yVs + (d - 3)-Vm. 

y 



(3.14) 



(3.15) 



(3.16) 



(3.17) 



(3.18) 



In order to solve the previous system, the derivatives of the different coefficients are required. 
To keep here the exposition simple, these expressions are given in Appendix C. 
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Decoupling of the perturbations. In this paragraph we show how the perturbations Mi 
and 5*1 can be decoupled from one another. Although outside the mail line of the paper, the 
analysis presented here may provide a different way to understand the perturbations to the 
critical solutions. In any case the paragraph can be skipped in a first reading. 

Going back to the equations for the perturbation (13.51) . we can formally integrate each of 
the equations to find the following system of integral equations for Mi and Si (keep in mind 
that the prime indicates derivation with respect to log z) 

flog z 



MUz) 



/logz 
dxBM{x)S 
-oo 



SAz) = e^— 



SA-oo) + 



logz 



dxAs{x)Mi{x) e-/-^'^-'«s(-0 



(3.19) 



Using again (13.51) they become 

/logz 
dx Bm{x)Si{x) e 
-oo 

flogz 



S[ - BsSi 



Si(-oo) + 



dxAs{x)Mi{x)e-r-^'^'^-^^'^ = e^-r'^-^sW 



— oo 



M[ - Am Ml 
Bm 



. (3.20) 



Taking now a derivative with respect to logz on both equations, and after some manipula- 
tions, we arrive at a set of two decoupled second-order linear differential equations for the 
perturbations Mi and 5*1 



M(' + VmM[ + Qm Ml 
S'l + VsS[ + QsSi 



0, 



0, 



(3.21) 



where the coefficients are given by 

Vm = - 

Vm = - 

Qm 



Am + Bs + ^ 
Am + Bs + ^ 



B[ 



M 



AmBs — AsBm — A'm + Am-;^ — , 



Qs = AmBs-AsBm-B's + Bs ' 



As' 



(3.22) 
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The system of equations fl3.2ip can be written in a more suggestive form by a further 
function redefinition 

Si{z) = 7/'5(z)e-5/-r'^-^s(x)^ (3_23) 
where the new functions iPm{z) and ipsi^) satisfy one-dimensional Schrodinger equations 

V'm + ^m^m = 0, 
^m + Vs^Ijs = 0. (3.24) 
The potentials for these equations are given in terms of the coefficients of Eq. (13.211) by 

Vm = Qm - ^V'm - -^M, 

Vm = Qs-\v's-\vl. (3.25) 

Hence, we can obtain the perturbations to the self-similar background by solving a system of 
two Schrodinger equations with the appropriate boundary conditions. The potential associated 
with these equations is completely determined by the self-similar solution. Once (13.241) is 
solved, the perturbations for the other functions, yi{z) and ^7i(-z), can be written in terms of 
iI)m{z) and ipsiz)- This formulation of the pertubations can provide a good arena for qualitative 
analysis. We will however proceed with a direct study of the original equations. 

Asymptotic behavior of the perturbations near z = and z = oo. In order to solve 
the perturbations Mi(z), Si{z) and yi{z) by integration of the system of differential equations 
(13.51) starting at z = we need to determine the behavior of these functions at the origin. To 
do so, we start by noticing that the function y(t,r) (including the perturbation) should have 
the limiting value determined in Eq. (12.301) . Because this condition is satisfied already by the 
self-similar background solution we conclude that the perturbation yi{z) has to vanish in this 
limit 

Z/i(0+)=0. (3.26) 
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Hence, studying the system fl3.2p around z = one arrives at the identity 

M[{z) + k{d-l)S[{z) = -A 

For positive A (as it is the case of the growing mode we are interested in), the integration 
of this equation resuhs in the combination Mi{z) + k{d — l)Si{z) blowing up at 2; ^ 0^ as 
2;"'*'. Therefore, in order to impose regularity of the perturbations at the origin we require the 
numerical coefficient of this divergent terms to be zero. This leads to the condition 

= 0. (3.28) 

Equations (13.26^ and (13.28^ give the initial conditions for the numerical integration of the 
system of differential equations (13. 2p at 2; = 0. Notice that since we have a linear system, only 
the ration Mi/ Si at the origin is relevant. Furthermore, the same conclusions can be reached 
from the equations determining yi from Mi and Si together with standard asymptotic analysis 
at 2; = 0. 

The functions Mi{z) and 5*1(2:) can be written, around z = 0, as a power series. From the 
reduced system of equations (I3.15P it can be seen that they can be expressed as 

00 

n=l 

00 

^W/^TH^rni+TT. (3.29) 

n=l 

Comparing with Eq. (12.321) we learn that the expansion variable here is the same as the one 
for the self-similar background solution. 

We now consider the behavior of the perturbations at infinity. Since in that limit the 
function y{z) tends to one, we find that the system of equations (13.21) simplifies for large z to 

{l + k)M[ = {3 - d)kyi - X[Mi + {d - 3)kSi], 

{l + k)S[ = yi-^[Mi + {d-3)kSi], (3.30) 

leading to 

Mi{z) + k{d-3)Si{z) 
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Mi{z) + k{d~l)Si{z) 



(.^0+). 



(3.27) 



Mi{z) + k{d-l)Si{z) 



-A 



Mi{z) + k{d-3)Si{z) 



(3.31) 



For the growing mode, A > 0, this means that the combination Mi{z) + k{d — 3)Si{z) vanishes 
at infinity as 



Mi{z) + k{d-3)Si{z) z-^ (z^oo). (3.32) 

To get the individual behavior of Mi (z) and 5*1 (z) we look at the asymptotic form of yi (z) 
for large z. From Eq. (13.31) . using the asymptotic expansion of the background functions when 
2; — > oo, we arrive at 

yi{z) ~ \Mi{z) -{d- 3)Si{z)\ (when z 00). (3.33) 

Plugging this in Eq. (13.21) we get the following equations, valid for large z, 

M[{z) = -AMi(z), 

S[{z) = -\Si{z), (3.34) 

implying 

Ml (2) = Mi'^^z'^ + ... 

Si{z) = S[°°^z-^ + ... (3.35) 

To summarize, we have found that for A > (the growing mode) all perturbations Mi{z), 
Si{z) and yi{z) go to zero at infinity as powers of z. On the other hand, for small z the 
only perturbation that vanishes is yi{z), while Mi{z) and Si{z) tend to constants in that limit 
satisfying the condition fl3.28p . 

The gauge mode. In studying the perturbations to the self-similar solutions of the Einstein 
equations it is very important to eliminate those perturbations which can be written as mere 
gauge transformations of the background solutions. In particular, by performing changes of 
coordinates on the self-similar solutions it is possible to introduce an explicit dependence on 
time that is just a gauge artifact. One can consider reparametrizations in time [lOj, leaving 
invariant the form of the line element (12. ip . Infinitesimally 

r — >r, t — >t + e{t). (3.36) 
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This changes the scahng coordinate z by 



r — > r + (5( 



r 



1 + 5(r) 



(3.37) 



where the function 5(r) is defined in terms of e(f) as 



(3.38) 



The functions M(t, r), S'(t, r) and ?/(/!:, r) only change under fl3.36p through their dependence 
on t. In the case of the self-similar solutions we find 



Z{z) — . Z{z) + b(T)Z\z) = Z{z)[l + 6{T)[logZ{z)]'], 



(3.39) 



where Z stands for M, S and y. This has precisely the form of the perturbation (13. ip if we 
choose the arbitrary function 6{t) of the form 

A 



5{T) = e-e'\ (3.40) 
with A any real constant. The perturbations Mi{z), Si{z) and yi{z) are then given by 

M^iz)„ = = log M{z) 

'S'l(2;)gause = ^ log S{z) 



jauge 



A 
2" 
A 

A 
2" 



(3.41) 



log?/(z) 



Using the behavior of the self-similar solutions near z = it is straightforward to see 
that the perturbations (I3.4ip automatically satisfy the boundary condition fl3.28p . However, 
in order to satisfy the boundary condition at z = cxd we find from (I2.33p that the Lyapunov 
exponent A has to be 

1-k 



Agauge- 

In particular, we see that the gauge mode is independent of the space-time dimension. 



(3.42) 
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3.2 Numerical study of the perturbations and calculation of the 
Choptuik exponent 

As in the case of the self-similar solution, the equations for the perturbations have to be 
numerically integrated using the background solutions obtained numerically in Sec. 12.31 We 
use a strategy similar to the one applied there. We numerically integrate the system (13. 5p 
starting at z = using a fourth-order Runge-Kutta routine. We find that for generic values 
of the Lyapunov exponent A the perturbation is singular both at the sonic point and at the 
value z = Za where y{za) = 1. We then scan in A to find a critical value for which the solution 
smoothly crosses both points. At the same time we ignore the value of A associated with the 
gauge mode. The singularities of the equations lead to some numerical instabilities, but the use 
of I'Hopital rule explained above gives a guidance to the numerical integration. A good deal 
of work leads to the results for the Choptuik exponents in Tabled] for dimensions d = 4,5,6 
and 7. As an important check of our results we agree with previous results in ci = 4 flU{ ITT] . 

In Figs. [7] and [8] we show the plots of Mi, Si and yi versus log 2; for three different values 
of k. As expected, all perturbations go to zero at large z, whereas only yi vanishes in the limit 
z —>■ 0^. This means that, for constant r, self-similarity is mostly broken near z = whereas 
the symmetry is restored at large values of z. In particular, we see that the position of the 
sonic point Zsp (marked by a dot in the plots in Figs. [7] and [H]) approximately indicates the 
value of z where Mi, Si and yi begin to decay towards zero. 

Notice, however, that the perturbations to the self-similar solution come multiplied by e"^'^ 
where A > is the positive Lyapunov exponent. This means that the region where CSS is 
broken grows as r increases until at some point the linear perturbation regime becomes invalid. 
The system goes then into a full nonlinear regime that, for supercritical solutions, results in 
the formation of a black hole in the asymptotic future. 

In Figs. [9] and [TO] we plot the perturbations Mi, Si, yi in terms of y. As in the case of the 
critical solution, this shows that the perturbations reach their asymptotic values at 2; = 00 
after a single turn. Once again this shows that we have correctly identified the most relevant 
perturbations leading to the computation of the critical exponents. 

The regularity of the perturbations at the sonic point fully determine the value of the 
Lyapunov exponent A. In Table [T] we show the values of both A and the Choptuik exponent 
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k 




^d=5 


^d=6 


^d=7 




k 


7(i=4 


7(i=5 


7d=6 


7d=7 


0.01 


8.747 


4.435 


3.453 


3.026 




0.01 


0.114 


0.225 


0.290 


0.330 


0.02 


8.140 


4.288 


3.376 


2.974 




0.02 


0.123 


0.233 


0.296 


0.336 


0.03 


7.617 


4.152 


3.302 


2.924 




0.03 


0.131 


0.241 


0.303 


0.342 


0.04 


7.163 


4.027 


3.233 


2.876 




0.04 


0.140 


0.248 


0.309 


0.348 


0.05 


6.764 


3.911 


3.169 


2.831 




0.05 


0.148 


0.256 


0.316 


0.353 


0.06 


6.412 


3.804 


3.107 


2.788 




0.06 


0.156 


0.263 


0.322 


0.359 


0.07 


6.099 


3.703 


3.049 


2.746 




0.07 


0.164 


0.270 


0.328 


0.364 


0.08 


5.818 


3.609 


2.993 


2.706 




0.08 


0.172 


0.277 


0.334 


0.369 


0.09 


5.565 


3.521 


2.940 


2.668 




0.09 


0.180 


0.284 


0.340 


0.375 


0.10 


5.334 


3.438 


2.890 


2.631 




0.10 


0.187 


0.291 


0.346 


0.380 


0.11 


5.124 


3.360 


2.841 


2.595 




0.11 


0.195 


0.298 


0.352 


0.385 


0.12 


4.932 


3.286 


2.795 


2.561 




0.12 


0.203 


0.304 


0.358 


0.390 


0.13 


4.756 


3.216 


2.751 


2.527 




0.13 


0.210 


0.311 


0.364 


0.396 


0.14 


4.593 


3.149 


2.708 


2.494 




0.14 


0.218 


0.318 


0.369 


0.401 


0.15 


4.442 


3.086 


2.667 


2.464 




0.15 


0.225 


0.324 


0.375 


0.406 


0.16 


4.301 


3.026 


2.627 


2.433 




0.16 


0.232 


0.330 


0.381 


0.411 


0.17 


4.170 


2.968 


2.589 


2.414 




0.17 


0.240 


0.337 


0.386 


0.416 


0.18 


4.048 


2.913 


2.552 


2.377 




0.18 


0.247 


0.343 


0.392 


0.421 


0.19 


3.933 


2.860 


2.517 


2.348 




0.19 


0.254 


0.347 


0.397 


0.426 


0.20 


3.825 


2.809 


2.482 


2.321 




0.20 


0.261 


0.356 


0.403 


0.431 


0.21 


3.723 


2.760 


2.449 


2.297 




0.21 


0.259 


0.362 


0.408 


0.435 


0.22 


3.627 


2.713 


2.417 


2.272 




0.22 


0.276 


0.368 


0.414 


0.440 


0.23 


3.536 


2.668 


2.386 


2.246 




0.23 


0.283 


0.375 


0.419 


0.445 


0.24 


3.449 


2.625 


2.355 


2.224 




0.24 


0.290 


0.381 


0.425 


0.450 


0.25 


3.367 


2.583 


2.325 


2.202 




0.25 


0.297 


0.387 


0.430 


0.454 



Table 1: Values of the Lyapunov exponent (left) and the Choptuik exponent (right) as a 
function of k for d =4, 5, 6 and 7. All numbers are calculated with a precision ±0.001. 
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7 = A^^ for various values of k in four, five, six and seven space-time dimensions. All values 
shown in this table are given with a precision of ±0.001. These values of the Lyapunov and 
Choptuik exponents are the main concrete result of this paper. Looking at the value obtained 
for the conformal fluid in five dimensions k = 1/4, X = 2.58, we find it is reasonably close to 
the value for a massless scalar field. This situation is very similar to the four- dimensional case, 
where the Choptuik exponent for a conformal fluid is very close to the value of a massless scalar 
field. This is also the case for the conformal fluids in c? = 6 (fc = 1/5) and d = 7 {k = 1/6), 
as compared to the values of the Choptuik exponents for the collapse of a massless scalar field 
in these dimensions [5]. Furthermore the value of the Choptuik exponent for the conformal 
perfect fluid in d = 5 is also close to the BFKL exponent as calculated in [Ij. To what extent 
this numerical agreement is serendipity or a sign of something deep will be analyzed in the 
next Section. 
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GRHA-dim = 5 
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Figure 7: Plot of the perturbations Mi{z) and Si{z) versus logz for k — 0.01, k — 0.1 and 
k — 0.2 in five dimensions. The dot gives the position of the sonic point. 
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GRHA-dim = 5 
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-10 -5 5 10 

Log z 

Figure 8: Plot of the perturbation yi{z) as a function of log 2;, both for k = 0.01, k — 0.1 and 
k = 0.2. As usual the dot indicates the position of the sonic point on each curve. 
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Figure 9: Plot of Mi versus yforA; = 0.01, A; = 0.1 and A; = 0.2 in five dimensions. 
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GRHA - dim = 5 
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4 A holographic glimpse through the Regge region 
4.1 Nonperturbative physics at weak 't Hooft coupling 

Probably the most important lesson we have learned from Maldacena's conjecture [19j is that 
the QCD string, i.e. the effective string defining the representation of Yang-Mills Wilson 
loops as sums over random surfaces, is the fundamental string in some higher dimensional 
geometry |20j. In spite of many efforts during the last ten years, we are still lacking the 
concrete string background geometry for QCD. However, many of the results so far obtained 
hint at the existence of a close relation between nonperturbative effects in QCD and black hole 
dynamics in the gravity dual picture. 

An interesting perturbative regime where stringy aspects of QCD naturally appear is the 
Regge limit of scattering amplitudes. In this region the center of mass energy is much larger 
than all the other Mandelstam invariants in the process under consideration, s ^ —t. There 
also exists some hard scale, ^ ^qcd' which allows for a perturbative treatment of the 
problem. Although the coupling is small, the amplitudes are dominated by logarithms of s, 
making a resummation of terms of the form [a^ log(s/(5^)]" mandatory. Note that ag = c^sNc/t^ 
is the usual 't Hooft coupling. This resummation can be performed using the BFKL ^14j 
formalism which predicts a Regge behavior of the type s"^^*-* for total cross sections, with uj{t) 
being some effective Regge trajectory. When the quantum numbers exchanged in the t-channel 
are those of the vacuum, the dominant trajectory is that of the BFKL Pomeron (also known 
as hard Pomeron since it is of perturbative origin). 

In the leading logarithmic approximation (LLA), where coupling and large logarithms carry 
the same power in the resummed terms, at asymptotically large energies and in the physical 
region of negative t the BFKL Pomeron trajectory goes like uj{t) ~ 0:^4 log 2. This power-like 
behavior of the cross sections calls for modification since it violates unitarity bounds at very 
high energies. In the context of deep inelastic scattering (DIS) and evolution equations for par- 
ton distribution functions, the onset of unitarity can be interpreted in terms of "saturation" 
corrections which tame the power-like growth with s. The physics of saturation is driven by 
perturbative degrees of freedom, because we are far away from the confinement region. How- 
ever, the treatment of the problem is nonperturbative since it corresponds to the description 
of a very dense system. The main target of this Section is to provide a holographic or gravity 
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dual description of these non perturbative saturation effects. 

Tlie gravity dual interpretation of the BFKL Pomeron has been previously addressed in [2] . 
Using a string target space in the ultraviolet of type AdSs it is possible to study the hard 
Pomeron at strong 't Hooft coupling cis ^ 1. This is done by computing perturbative string 
amplitudes in the Regge limit with log s being the small perturbation parameter, where 
is the string scale, R the curvature radius and C-H = ^ 1. In this weak gravity limit 

the BFKL diffusion in transverse momentum is mapped into propagation along the curved 
holographic direction. The physics of diffusion in logs time at weak and strong 't Hooft 
coupling is qualitatively the same but with different coefficients. 

Let us remark that in this holographic correspondence the BFKL Green's function in the 
weak gravity side, computed in a small curvature gravitational background, is mapped to the 
large as BFKL Green's function in the gauge theory side computed around the saddle point 
"background" 7 = 1/2 of the BFKL kernel eigenvalue x(7) = 2\I^(1) - \I^(7) - \E'(1 - 7). As 
we will explain in the next Subsections, the nonperturbative phenomena of saturation requires 
to move to a different saddle point background 7^ > 1/2 where amplitudes are not power- 
like anymore at large s and enjoy "geometrical scaling". In the gravity dual picture this 
change of saddle point "BFKL" background should correspond to a nonperturbative change 
of gravitational background. This new background would require to take into account off shell 
gravity effects and therefore is beyond the semi classical weak gravity limit. 

In Ref. [21] multi-Regge kinematics was extended to include nonperturbative effects in 
the region of small momentum transfers of the BFKL equation in a c? = 4 string theory. In 
Ref. [22] , assuming that the hard Pomeron intercept in A/" = 4 SYM equals the graviton one in 
AdSs X S^, it was possible to obtain a very good extrapolation of DIS anomalous dimensions 
of twist-2 operators from weak to strong coupling. 

After this brief introduction we now move to a more detailed description of our holographic 
picture. We start with an introduction to the calculation of scattering amplitudes in the Regge 
limit of perturbative QCD. 

4.2 Deep inelastic scattering and BFKL 

Let us start with a short introduction to DIS at small values of Bjorken x. In very simple terms, 
in DIS we are resolving the structure of a large target, e.g. a hadron, with an electromagnetic 
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probe. This probe is typically a photon with virtuality and transverse size In the 

target's rest frame we can think of the photon's wave function as forming perturbatively well 
before the interaction with the target. When the transverse size of the photon is small it goes 
through the target like a "needle", however, for smaller virtualities the leading component of 
the photon's wave function, a quark-antiquark pair, is dressed with a color cloud of higher 
order gluonic components. In this case the 7*-hadron cross section behaves very similarly to 
that of hadron-hadron scattering (see Ref . [23] for a review) . 

In a leading twist picture, the virtual photon resolves partons in the hadron carrying a 
fraction x of its longitudinal momentum. This takes place with a probability g{x,Q'^), the 
so-called parton distribution function. The evolution with of g{x, Q^) is well described in 
perturbation theory using the DGLAP formalism [23]. As the hadron is a complicated strong 
interacting object it is not possible to calculate the parton distribution functions from first 
principles. Nevertheless, very good fits to experimental data have been achieved by assuming 
a certain nonperturbative functional form for the x dependence of the parton distribution at a 

close to the confinement scale and evolving it to larger scales using perturbative evolution 
equations. CoUinear factorization theorems ensure that, if the final is large enough, these 
parton distribution functions are universal: they can be measured in DIS experiments in 
electron-hadron machines and then used for predictions at hadron-hadron colliders. The region 
of a; 1 is well understood from the experimental and theoretical sides. 

The situation in the limit a; — > is far more complicated. In principle, even if is large 
enough as to ensure a small value for the strong coupling, it turns out that the cross sections 
or structure functions are dominated by logarithms of x which are very large as we enter this 
region. This implies that terms of the form log (l/x) should be resummed to all orders to get 
an accurate description of the data. The resummation of these terms can be effectively carried 
out by means of the BFKL [14j evolution equation, which can be considered as orthogonal 
to DGLAP in the (x, Q"^) kinematic plane. While DGLAP is based on perturbative splittings 
where one emission is ordered in transverse momentum with respect to the previous one, BFKL 
stems from strong ordering in longitudinal components. After integration over the phase space 
of these emissions a power like growth of the form x~^ is obtained for g(x,Q^). When the 
resummation is of leading logarithms in x, i. e. if the power in the coupling is the same as that 
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of the logarithm then, for x ^ 



A 



asN, 



c 



4 log 2. 



(4.1) 



TT 



When one power of the logarithm is lost we are in the next-to-leading order approximation 
and A is smaller. 

The BFKL equation in LLA has been rederived in many different ways in the literature. 
The original calculation was based upon the study of scattering amplitudes in the so-called 
multi-Regge kinematics. There are two main ingredients: the elastic scattering amplitude for 
the process g + g g + g and the inelastic amplitude for g + g—*g + g + g. We focus on 
gluons since their splitting function dominates at small x. 



Figure 11: Four point amplitude represented as the exchange of a Reggeized gluon in the Regge 
limit s ^ —t. 

When we consider the M2^2 scattering amplitude in the Regge limit of s ^ — t it has a 
simple factorized form in terms of two external t-dependent functions and an energy dependence 



This exponentiation of the logarithms in energy is known as Reggeization of the gluon (see 
Fig. [TTl) . If we project on the gluon's quantum numbers in the t-channel this means that, while 
at lowest order in as we have the usual gluon propagator, when higher orders are taken into 
account this propagator is modified by the Regge factor above-mentioned. The gluon then lies 
on a Regge trajectory, j(t) = 1 + ujQ^t), which, as usual, reproduces the gluon's spin when t is 
equal to the gluon mass, i.e. j{t = 0) = 1. In more detail, we have 



t 



s 




of the form (^)'^°^*\ In DIS, X ~ Q"^/ s and the limit of large s corresponds to small x physics. 
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where t = —q 



^ and n is an infrared regulator. 




Figure 12: The BFKL Pomeron constructed from two Reggeized gluons exchanged in the t- 
channeL Without interactions between them this amphtude is infrared divergent and decreases 
very rapidly with s. q is the momentum transfer. 

It is possible to construct a perturbative Pomeron with the quantum numbers of the vacuum 
if we exchange two Reggeized gluons and project the amplitude on a color singlet. In this case 
the four-point amplitude can be written as 



where fcj=A,B are the two-dimensional transverse momenta of the incoming gluons. Note that 
the transverse components of the four momenta decouple from the longitudinal ones. These 
are integrated out in multi-Regge kinematics with the Bjorken variable x playing the role of a 
cutoff. Thus we can consider Y ~ logs ~ log as the cutoff of the effective theory for the 
transversal degrees of freedom and the BFKL equation as the renormalization group evolution 
in "logs" time. 

Let us indicate that Eq. (14.31) corresponds to the amputated gluon Green's function since 
for zero Y we get the normalization 



which corresponds to the t-channel exchange of two normal gluons. 

It is interesting to remark some simple facts derived from Eq. (14.31) . The trajectory in 
Eq. (14.21) is a negative function, which means that the elastic scattering amplitude in Eq. (14. 3 p 
decreases as the energy grows. The delta function in Eq. (14. 3 p imposes conservation of trans- 




(4.3) 




(4.4) 
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verse momenta but it also means that we can consider the evolution of a Reggeized gluon as 
the action of a diagonal operator acting on some basis in the transverse momentum space. 

Being a color singlet, the BFKL Pomeron does not radiate and it is responsible for "diffrac- 
tive" events, where there are regions in the detector without hadronic activity, the so-called 
"rapidity gaps". To make connection with total cross sections or structure functions we can 
use the optical theorem. In this way we can relate the imaginary part of the forward elastic 
amplitude with the total cross section. This means that we can use Eq. (14.31) with g = to 
obtain a first contribution to g{x, Q"^). 




Figure 13: Ladder-like structure for the BFKL Pomeron. The interactions between the two 
Reggeized gluons are built with gauge invariant effective vertices. The amplitude is now in- 
frared safe. 



As it stands, Eq. (14.31) is not an infrared safe quantity since it has an explicit dependence on 
/i. This dependence is eliminated when we include the second ingredient in the BFKL equation: 
the inelastic A^2^3 amplitude. This amplitude is calculated in multi-Regge kinematics where 
the three outgoing gluons are well separated in rapidity from each other. This gives rise to a 
gauge invariant effective emission vertex. When this vertex is included we have that the gluon 
Green's function in the forward case now reads 



e2^o(-fci)Y 



5(2) ( k 



kf 



n=l i=l 



'^e (kl - 1% e'4^'-''y^s^^) [k^-k^ + Y^k, 



where 



kA + ^ki 



1=1 




(4.5) 



(4.6) 
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and yo = Y (see Fig. [T3|) . Below we will show a simpler representation where the dependence 
on the infrared regulator explicitly disappears. Before this we would like to point out that the 
action of the real emission kernel is not diagonal in transverse space anymore, as can be clearly 
seen in the second line of Eq. (14. 5p . This will be an important point when connecting with 
gravitational theories. It is also worth noting that the Green's function without virtual terms 
would be 

+ E ^ n / (e - ^0 {'-^ - + E ' 

n=l ' i=l ^^i \ l=\ / 

this means that the real emission makes g{x^ Q"^) grow very rapidly as we move towards smaller 
values of x. Let us note that both competing effects in Eq. (14.31) and Eq. (14. 7p depend on the 
infrared regulator /i. 
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Figure 14: Dependence on 7 of the eigenvalue of the BFKL kernel for different values of n. 

Making use of bootstrap, the virtual and real contributions to the BFKL amplitude can be 
combined in such a way that the infrared parameter /i is eliminated. This brings the following 
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simpler representation for the Green's function 



J 2m \kl I cu - asXnil) ' 



where 



X„(7) = 2vl/(l)-vl>(^^ + ^j _vl/(^l_^ + ^j (4.9) 

is the eigenvalue of the BFKL kernel and ^ is the logarithmic derivative of Euler's Gamma 
function. Q is the azimuthal angle formed by the two-dimensional vectors kA and ks- Its 
Fourier conjugate variable n can be interpreted as a conformal spin in elastic scattering. The 
dominant region of integration is 7 ~ 1/2. Around this point we can see in Fig. [TUthat the 
eigenvalue of the kernel is positive only for n = 0. This implies that the dependence on n is 
only important at low energies and we can keep n = for our arguments. Hence, we average 
over the azimuthal angle 

(^(fcA,fci3,Y) = (f(^kA,kB,Y^ 

1 fdufdj fk\y-'^ 



T^kAks J 27ri J 2ni Xk"^ J u — agxil)^ 

where x = Xo and, from now on, ki = \ki\. 

Let us now study the structure of this equation for very large values of a^Y. We can write 

'fik..ks.Y) = -A- f #y'^*"-^ (4.11) 



X(',) ^ 41og2+ 14(3(7-^1 +■•• (4.12) 



TikAks J 2Tii \k\ ^ 

and evaluate the integral using the saddle point approximation around the minimum of x at 

7 = 1/2 

to obtain the expression 

i^ikA.kB^Y) ^ -^^e^^-=4==e^^, (4.13) 
2-nkAkB Vl47rC3asY 

with t = \og{k\/k\). This asymptotic expansion is driven by the first exponential, which 
contains the leading behavior in energy, with A being the Pomeron intercept given in Eq. fl4.ll) . 
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The function xil) has poles at integer values of 7 with the physical region corresponding to 
< 7 < 1. The poles at 7 ~ 0, 1 are related to virtualities in the internal propagators smaller 
and larger, respectively, than the external scales. The terms containing (3 are sensitive to 
the shape of x its minimum and indicate the speed of propagation into the infrared and 
ultraviolet regions. It is easy to verify that the function 

^{kA,kB,Y) = kAkB^ikA,kB,Y) (4.14) 

asymptotically fulfills the diffusion equation 

4 1og2$ + 14C3i^. (4.15) 



This diffusion in transverse momentum is IR/UV symmetric as a consequence of the invariance 
of X under the transformation 7 — 1 — 7. 

The gauge/string duality has offered the possibility to explore the ^ 1 sector of DIS 
with the hope of building the correct string dual of QCD. In the dipole approximation to DIS 
a holographic description at strong 't Hooft coupling was obtained in terms of a semiclassical 
approximation to the light-like Wilson loop in the target metric. At finite temperature and 
in the quark-gluon phase this was done using the AdS black hole metric in [25j. Another 
interesting study was carried out in Ref. [26j where structure functions were calculated at 
strong 't Hooft coupling. The model is based on A/" = 4 supersymmetric Yang-Mills (SYM) 
theory broken down to A/" = 1 at a given mass scale. When going to the small x limit the 
analysis was performed within the single string picture. This set up is limited since, as we 
will see in the coming Subsections, new effects in the gauge theory side related to the onset 
of unitarity indicate that one should better operate with multi-string configurations. A more 
complete analysis of this region was presented in Ref. [3|. In that work, rather than resumming 
multi-loop string amplitudes, inspiration is taken from perturbative QCD to understand the 
flow towards unitarity at large as and the gravity effect on a "bulk" photon of large virtuality 
is worked out semiclassically. This photon is defined using any of the U{1) gauge fields in the 
bulk associated with the global /^-symmetries of the theory. With increasing virtualities of this 
probe the transition to unitarity changes from a single-Pomeron picture to one with multiple 
graviton exchanges. 



54 



In the present work we investigate unitarization effects in DIS in tlie limit of weak 't Hooft 
coupling (is and very large center of mass energy. In this limit the LLA, where terms of the 
form (agY)"' are resummed, is accurate since higher order corrections are always suppressed by 
larger powers in the coupling, e.g. the next-to-leading order approximation (NLLA) resums 
as{asY)"' terms. The NLLA corrections in A/" = 4 SYM have been studied in [27]. The program 
to construct even higher order corrections in this theory was initiated in |28j . 

It is important to note that we are not working with a supersymmetric theory from the 
beginning. It is the dynamics of our process, treated in multi-Regge kinematics, what provides 
a highly symmetric effective theory at high energies. As a matter of fact, in the LLA, the 
BFKL kernel is the same in QCD or A/" = 1,2,4 super Yang-Mills theories. This is because 
fermions are suppressed at high energies and their effects are only important at higher orders 
in the approximation. Moreover, we are not sensitive to the running of the coupling since the 
Feynman diagrams contributing to it do not belong to the LLA because they bring an extra 
power in as without generating a logarithm in energy. Last, but not least, this theory naturally 
shows planarity which comes as a consequence of the LLA. 

We follow a bottom-up approach and find that the perturbative description of the satura- 
tion line, to be defined in detail below, within the BFKL formalism shares many interesting 
features with the critical formation of a tiny black hole in higher dimensions. Our approach 
to the identification of the Regge holography will be mostly constructive. The first issue is to 
identify the holographic extra dimension. Since in the BFKL approach to the Regge limit we 
are working with an effective theory parametrized by logs, we define the holographic coordi- 
nate as a^Y ~ (is logs ~ Og log the natural energy scale in the problem. In the mapping 
of BFKL dynamics to critical black hole formation we consider the evolution in the holographic 
variable a^Y as dual to time evolution in the gravitational collapse. At the threshold of for- 
mation of the black hole the variations in the metric components are large and should be well 
described within a perturbative analysis in the gauge theory side. 

The next step in our description should be to identify the geometry on the holographic 
direction. This geometrical information for Regge holography should be hidden in the BFKL 
kernel x{l)y with 7 parameterizing a particular gravitational background. In the AdS/CFT 
correspondence, in order to describe a phase transition, we need to deal with the full quantum 
gravity canonical ensemble of metrics with a given boundary topology. In our case, to describe 
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the transit through the saturation region, which corresponds to the flow from a dilute partonic 
regime to a dense system, we need to consider the full ensemble of BFKL backgrounds weighted 
by the kernel xil)- the next Subsection we will try to put some substance to these general 
guiding lines of work. 

4.3 Unitarity in DIS and perturbative saturation 

It is well known that the LLA in Eq. fl4.13p . predicting a power-like growth in energy for total 
cross sections, violates the unitarity bound (logs)^ set by Froissart. The problem of unitariza- 
tion of cross sections at high energies is a very complicated one. From a perturbative point of 
view it can be understood by introducing a new set of Feynman diagrams. In the context of 
QCD Reggeon fleld theory a new vertex arises: that of the transition from 2—^4 Reggeized glu- 
ons. This "triple Pomeron vertex" is invariant under SL{2,C) transformations [2S], a Mobius 
invariance also present in the BFKL Hamiltonian in the non- forward case [30] . In the large Nc 
limit the triple Pomeron vertex in QCD reduces to the vertex used in the Balitsky-Kovchegov 
(BK) equation [TB]. This equation describes a DIS process where the target is an extended 
and dense object. The effect of taking into account unitarity in this equation is to introduce 
a non-linear piece added to the usual BFKL kernel. Neglecting impact parameter effects and 
azimuthal angle dependence, the BK equation can be written as 

d<l>{kA,kB,Y) 



^'^ dx 



+ 







1 — X 

with the initial condition flxed by 



= -^{kA,kB,Yy 

$ {V^kA, kB, Y) + 1$ (^-^, ^B, - 2$ {kA, kB, Y) 



(4.16) 



1 fdr^f}^'^'- 
7T / 2m \ kl 



$(/cA,fc5,Y = 0) = - -L-A-^] . (4.17) 



The non-linear term in Eq. (14.161) arises from the introduction of fan diagrams and it has a 
mild effect at the beginning of the evolution, when the gluon amplitude is small. The growth is 
therefore power-like at early times in a^Y. On the contrary, for large values of a^Y, when the 
amplitude is larger, unitarity appears in the form of nonlinearities and the quadratic term in 
Eq. (I4.16P becomes very important since it completely suppresses the growth of the amplitude 
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with energy. In terms of parton densities this transition to a saturated state can be understood 
as the crossover from a dilute system of partons within the hadron to a very dense one which 
can be interpreted as a "color glass condensate" (CGC) [31] . 

As we already remarked, the precise dynamics of unitarization requires to take into account 
multi-Pomeron vertex effects. However, it is still possible to describe within a single Pomeron 
picture some general features such as scaling laws present in the crossover to the saturation 
region. For this we closely follow [32] . 

We are considering the limit of very large energies and we can start by evaluating the integral 
in Eq. (14. lip asymptotically. It is important to correctly identify the virtuality corresponding 
to the target in the DIS process, k\ = Qtarg Eq. (14.111) . and that of the projectile, /c^ = Qproj- 
We can then write 



^(Qtarg,Qproj,Y) = -7^/^ ^ e'^^^^-^ (4.18) 

At a^Y 3> 1 it is possible, for fixed values of Qproj = Qo ^^"^ ^^^^ 7-background 7 = 
7o(Qo5 Y) that dominates asymptotically. It is implicitly given by the saddle point condition 

x'(7o)«.Y + log(^) = 0. (4.19) 

To evaluate the Green's function we now expand around this saddle point 

X(7) ^ x(7o) + x'(7o)(7 - 7o) + ^x"(7o)(7 - 7o)' + • • • (4.20) 

and obtain 



^0 



p 2x"(7o)"sY 

TrQtargVX (7o)2vra,Y 

where to = log (Qo/'5proj)- The leading behavior of this expression is driven by the first 
exponential. The remaining pieces are of diffusive origin and do not need to be considered in a 
first analysis. The growth with energy of Eq. (I4.2ip is dictated by the function x{l) ~ 7X'(7) 
which we plot in Fig. [151 This function is monotonically decreasing as 7 grows. It has a zero 
at the point 

7cr = ^ 0.6275... (4.22) 

X l7crj 
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Figure 15: The function of 7 dominating the energy dependence of the gluon Green's function. 



For any 70 in the physical interval (0,1) the kernel x(7o) is positive. For 70 < 1/2 < 7cr we 
have that x'ilo) is negative and Eq. fl4.2ip rises very rapidly as Y gets larger. In the range 
1/2 < 7o < 7cr, x'ilo) changes sign but its value is not large enough as to stop the growth 
of the Green's function. However, for backgrounds with 70 > 7cr there is no growth because 
x'ilo) is so large that it forces the amplitude to decrease with energy. The value at which this 
change in behavior takes place is x'ilcr) — 4.8833... 

The value 70 = 1/2 is precisely where the derivative of the kernel cancels, x'(l/2) = 0, 
and where the usual BFKL expansion happens. It generates a rise of the cross sections which 
must be tamed at some point in energy in order to respect unitarity bounds. Unitarity can be 
enforced in our single Pomeron picture if we impose the extra constraint that the dominant 
region of integration at asymptotic Y corresponds to 70 = 7cr- What we obtain is 



V-' iQtaigj Qproj! Y) 



27c: 



e2x"(7cr)<5sY 



Qproj / vrQt\,gA/x"(7cr)27ra,Y' 



(4.23) 



with tcr = log [<5cr(Y)/<5proj]; xilcr) — 3.0645... and x"ilcv) — 48.5176... In the kinematic plane 
of DIS, the critical line where the transition to saturation occurs can be read off the saddle 
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point condition in Eq. fl4.19p 



Qcr ( Y) — <5targ G^p 



X'(7cr) 



(4.24) 



This "saturation line" is represented in Fig. [161 Taking into account Eq. fl4.24p it is easy to 
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Figure 16: Diagram in the kinematic plane of DIS for the transition from a dilute partonic 
system to a dense one through BFKL evolution. The transition to a saturated phase takes 
place crossing a region with scale invariance. 



show that the saturated amplitude in Eq. (14.231) manifests an extra approximate symmetry 
at 7o = 7cr- The appearance of scaling symmetries is common in critical phenomena. In our 
case the BFKL Green's function with saturation is, up to very mild subleading corrections, 
invariant under the scaling transformation 



Qpro] 



asY + log A, 

Qproj 



targ 



(4.25) 



We can then identify x'(7cr)/2 — 2.4417... with a critical exponent. 

The invariance of scattering amplitudes under this scaling transformation is known as 
"geometric scaling" in DIS at small x and has been experimentally observed in HERA data 
for the 7*p cross sections in the region x < 0.01 over a large range in Q^, see Fig. [T71 
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Figure 17: HERA data for a^*p with x < 0.01 versus the variable r ~ Q'^x^^. 

4A Holographic map to critical gravitational collapse 

Now that we have described an approach to unitarization in high energy QCD we have all the 
ingredients needed to propose a holographic map with the formation of a tiny black hole. We 
identify the region of linear growth with a^Y in the QCD side with that of time evolution on a 
supercritical line just above the critical manifold of codimension one present in our gravitational 
studies of a perfect fluid (see Fig. [18]). We set the initial critical density distribution of perfect 
fluid at time t = — oo and evolve the system following Einstein's and matter equations of 
motion. 

It is also natural to map the radial component r in the gravity side to 1/Q in the gauge 
theory sector. Both variables measure the spread of the system on a space "orthogonal" to the 
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Figure 18: Trajectory in phase space of solutions for supercritical initial conditions. The 
transition to form the black hole is dominated by continuous self-similarity. 

holographic direction. Note that the typical virtuality should be large enough to ensure a 
perturbative description of the scattering. This brings us to the r ^ 1 region of the critical 
collapse, precisely the zone where the black hole is formed. 

The competition in QCD between Reggeon exchange, which, as we have already explained, 
tends to reduce the size of the amplitude, and the emission of on-shell gluons, with the opposite 
effect, has an interpretation during the collapse: it corresponds to the balance between self- 
gravitational attraction and kinetic energy, or pressure, in the perfect fluid. To understand this 
correspondence in more detail we can go back to Eq. fl4.22p and note that in the unitarized 
solution the dominant background is 7cr > 1/2. This is a non-trivial statement since now 
we are not at the minimum of xil) the original BFKL case. The dominant region of 

integration is shifted towards the ultraviolet pole at 7 = 1. This implies that the net effect of 
non-linear unitarity corrections is to suppress the diffusion towards the infrared region. 

This can be mapped into the dual gravitational collapse framework. The infrared modes in 
QCD have small virtualities and therefore occupy large regions in the transverse plane. The 
physics of saturation is that of a system with perturbative degrees of freedom forced to grow 
in a finite region. The critical line Qcr(Y) in Eq. fl4.24l) shows that, in order to accommodate 
further multiplicity, the preferred configurations are those with small transverse size, i.e. large 
virtualities. This is exactly what happens when the black hole is formed: the largest part of 
the initial collapsing matter does not confine in a compact size and expands infinitely, while 
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only a minute fraction remains within the tiny black hole radius in the region r < tbh- It is 
natural to think that the black hole contains the bulk of collapsing modes with the smallest 
spread in the radial component. 

At asymptotic energies we can heuristically picture the mechanism of suppression of infrared 
modes as that of a two-dimensional disk packing. It is useful to go back to the linear BFKL 
equation and notice that the evolution it describes is not sensitive to the transverse size of the 
gluons in the ladder. This is a consequence of the lack of ordering in transverse scales, which 
generates an IR/UV symmetric random walk in these variables. In some sense, the system does 
not have memory of the transverse size of the previous steps in the evolution and is insensitive 
to the size of the target. Note that 7 = 1/2, corresponding to the minimum of xil)^ is the 
only 7-background where the saddle point condition in Eq. fl4.19p does not mix longitudinal 
and transverse components. As it was shown in Eq. (14. 5p the gluon Green's function is the 
sum of many configurations with different multiplicities. We can then represent the asymptotic 



n) 



Figure 19: The evolution with energy of the linear BFKL equation is insensitive to the target's 
size and circles are allowed to overlap. 



linear evolution as a function of the multiplicity (n) in Fig. [19] where the average transverse 
size of each gluon is represented by the circles, which are allowed to overlap. The largest circle 
corresponds to the target's transverse size. 

Nearby the saturation line the parton evolution is constrained by the relation in Eq. fl4.24p 
which we can write in terms of transverse sizes T = 1/Q as 



T 



X'(7c 



-a.,Y 



(4.26) 



This implies that each gluon decreases its transverse size with respect to the previous one as 
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12 3 4 

Figure 20: When saturation is introduced the evolution with energy of the BFKL equation is 
sensitive to the target's size and circles are not allowed to overlap. 

the energy grows. This introduces ordering in kt and suppresses the random walks towards 
lower virtualities enhancing the diffusion into ultraviolet modes. The system can now resolve 
the transverse size of the target and it can be pictured as a problem of packing many hard 
disks within a circular boundary without allowing overlapping. This is pictorially shown in 
Fig. [201 The fact that we have CSS indicates that there is a fractal structure in the evolution 
since the system reproduces itself at shorter distances. 

In terms of the function in Fig. [15] we can associate gravitational potential energy with the 
term — 7x'(7) and kinetic energy with xil)- Thus, for 70 < 7cr kinetic energy dominates, with 
no formation of a black hole. This corresponds to violation of unitarity in the Yang-Mills side. 
The threshold for formation of a black hole corresponds to 70 = 7cr, when the gravitational 
potential energy is large enough to neutralize its kinetic counterpart. This is interpreted as 
the existence of a unitarization mechanism in QCD. 

Let us indicate again that the QCD transition from linear evolution, which corresponds to 
a "dilute" region in the transverse space, towards the "dense" regime is characterized by the 
critical anomalous dimension 7cr and the "saturation exponent", Acr = x'(7cr)/2 — 2.4416... It 
is remarkable that these are pure numbers independent of any kinematic variable or the 't Hooft 
coupling. The interesting phenomena discovered by Choptuik, which we have reproduced here 
for arbitrary dimension, is that the perfect fluid goes through a similar process at the threshold 
of formation of the black hole with its radius scaling as 

TBH ~ 4 b -Pi'/'"", (4.27) 
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with io some fiducial length scale. The exponent Abh is also a pure number which depends 
on the dimension and the speed of sound in the fluid, y/k. It is difficult to find compelling 
holographic arguments to select a particular value of k. Nevertheless, we believe that the limit 
of "radiation fluid" is the correct one. In this case k = l/{d — 1) and the energy-momentum 
tensor is traceless. This would support our intuition that a conformal symmetry underlies the 
crossover region towards the black hole formation. We find further evidence in this direction 
when we look at Tableland see that the values for A in the radiation fluid case are 2.58, 2.48 
and 2.42 in dimensions five, six and seven, respectively. These values are very similar to Acr in 
the gauge theory side. 

At the stage we are in our holographic analysis it is not possible to precisely determine 
in which dimension the dual gravitational physics must live. From the previous arguments it 
seems it should not be far from dimension four. In a previous analysis [Ij it was argued that 
the mapping could be to dimension five. This was based on the study of spherically symmetric 
collapse of a massless scalar field [7]. There, a relation identical to f l4.27p was found for the 
size of the tiny black hole formed when imploding an initial radial density of static scalar 
field. It turns out that in this case Abh is very similar to Acr in QCD when the scalar field 
propagates in five dimensions. As discussed in the Introduction, the collapse of a scalar field 
also manifests discrete self-similarity. This means that the metric and field components at 
the critical manifold, generally denoted by Z^{t, r), reproduce themselves after a finite echoing 
period A. This means that Z^{t,r) = Z^{e^t,e^r), where in five dimensions A 3.44 [3]. 
However, it turns out that there is no analogue of this echoing period in the Yang-Mills theory 
side. 

As we already have discussed, in the case of a spherical collapse of a perfect fluid with 
equation of state p = kp, there is also self-similarity but in this case it is continuous (CSS). 
This implies that we can write the critical solution as a function of the variable z defined in 
Eq. fl27Hl) only, i.e. Z,{t,r) = Z,{z). 

In Section [2] we found the critical CSS solution by solving the Einstein equations sourced 
by a perfect fluid. As an example in four dimensions, we can look at the profile of the function 
y{r,t) which is proportional to the ratio of the mean density inside the sphere of radius r 
to the local density at r. In Fig. HT] we show how y{r,t) maintains a constant r-profile for 
different values of t. This implies that the solution is CSS since any continuous change in 
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the time coordinate can be compensated by a change in r leaving the solution invariant. The 
mapping of this scaling symmetry in the QCD side corresponds to the dual geometrical scaling 
in Eq. fOSD . 

It is important to note that in QCD geometric scaling is only present in the region very 
close to the saturation line. In the same way the CSS present in the gravitational critical 
solution breaks down very rapidly when we introduce a Liapunov perturbation of the form 
shown in Eq. (13.11) . 

The Choptuik exponent characterizing the black hole radius of Eq. 04.271] can be obtained 
by searching for these Liapunov modes of instability of the CSS solution. In Fig. [22] we show 
how the exponentially growing mode removes CSS from the solution to the collapse. The rate 
of growth of this mode is given by a coefficient A which coincides with the Choptuik exponent. 
As explained in previous Sections, we have numerically extracted this coefficient in different 
dimensions. 

To summarize this Section, in QCD there exist two interesting types of nonperturbative 
physics: one is the confinement / deconfinement phase transition with a gravity dual [35] given 
by a Hawking-Page type of phase transitions [3l]. The other one, stressed in this article, 
is the phenomena of parton saturation or dilute/dense system of partons transition which 
we suggest to have its gravity dual in the critical formation of tiny black holes. Whereas 
in the confinement /deconfinement phase transition in Hawking-Page representation we have 
an entropic balance to characterize black hole nucleation, in the parton saturation case we 
have an infrared-ultraviolet energy competition which we map into a potential/kinetic energy 
competition in the gravity dual. 



65 





t = -0.10 


1 


1 1 1 




t = -0.12 








t = -0.14 


















///' 





-25 -20 -15 -10 -5 5 

Log (r) 

Figure 21: A typical solution to gravitational collapse of a perfect fluid with continuous self- 
similarity. 
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5 Unanswered questions and work in progress 



We are aware of the fact that we are exploring unchartered territory. On the gravitational side 
we are looking at the formation of small black holes, with the corresponding high curvatures 
and hence the need to introduce corrections in oi in the Einstein equations of motion. We try 
to relate holographically this region with a yet not fully understood region in QCD or Yang- 
Mills theory: the realm of the hard Pomeron. Some rudiments of the would-be correspondence 
have been spelt out in the previous Section. We list here a number of open questions, general 
philosophy, and work in progress. 

So far we have considered the collapse of either scalar fields or perfect fluids. The next step 
is to consider gravitational collapse in the case of a full type-IIB string theory. We are currently 
working on this problem, and there are some interesting provisos worth discussing. First the 
beta-function equations in type-IIB receive corrections starting at four-loop order in a', hence 
if we do not consider curvatures close to the string length our arguments should be reliable 
at least qualitatively. Within the spirit of Choptuik computations, since we are considering 
the formation of very small black holes, the cosmological constant should not be relevant in 
the computation of the Choptuik exponents. What is not clear however, is whether we are 
forming small five-dimesional AdS black holes J34j or a small ten-dimensional Schwarzschild 
black holes embedded in AdSs possibility shown in [36]. The Choptuik exponent is 

likely to depend on the various scenarios. This is currently being investigated, and we hope to 
report our results in the future [37j . 

In all our computations, supersymmetry does not seem to play any role. It is important 
to stress that the viability of the proposed holographic relation lies in the fact that we are 
computing critical exponents, normally protected by some form of universality. This is the 
only stability protection we can advocate against all corrections imaginable. In addition, on 
the gauge theory side we have an effective large- expansion generated by the BFKL analysis. 
Their effective A'efr is proportional to log s. To what extent this leads to a systematic "large- 
Nc expansion remains to be seen, but should such an expansion exist it would also provide 
some substance to our musings. 

The region we are trying to explore, at least on the gauge theory side is accessible to 
experiment. How the BFKL evolution is unitarized remains to be understood. Our claim is 
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that under a holographic map this should give precious clues on how gravitational collapse 
becomes "unitarized" before forming a singularity. On the QCD side the BFKL evolution is 
thought to lead to the color glass condensate [31] (for a review see [38] ) or some similar scenario 
in which we go from a single reggeized particle to a many-body state. This mechanism is 
somewhat reminiscent of the Klein paradox. 

In string theory [39] , it is possible to give a heuristic description of how, as we approach 
the formation of the horizon in gravitational collapse, the description in terms of a "single" 
string should evolve into a multi-string system, thus diluting the formation of a singularity. 
This analysis is based on entropy arguments, and for the time being we do not know how to 
holographically match these arguments with an entropic argument in the CGC The formation 
of the latter, as well as the scaling arguments in gravitational collapse are deeply dynamical 
phenomena, far from any notion of equilibrium. At this moment it does not seem easy to 
make a connection between both arguments. We believe, however, that the formation of the 
CGC (or whatever unitarizes the BFKL evolution) will, through a holographic map, provide 
a description of what happens in the final stages of gravitational collapse. If this picture is 
right, the exploration of the Regge region in future colliders is likely to give us a completely 
new panorama of what happens when black holes form. 
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Appendix A. The Einstein equations 



In this Appendix we detail the derivation of the Einstein equations for the hne element 



a{t, rydr + a{t, rydr^ + R{t, rydVLi 



(A.l] 



coupled to a perfect-fiuid energy-momentum tensor (12. 2p . Here dQ^^ is the line element for 
a (Z) — 2)-dimensional space of constant curvature k = 0, ±1. The case with k = 1 corresponds 
to the spherically symmetric case studied in this paper. 

In order to write the geometric side of the Einstein equations it is very convenient to 
introduce a system of tetrads 



3O = a{t, r)dt, = a{t, r)dr, e' = R{t, r)e' 



(A.2) 



where 6** is a system of tetrads for the transverse space satisfying dflf, = dijO'^O^ . 

The Riemannnian spin connection can be computed from the torsion-free conditions de"" + 



Lj"'i^ A = 0, with the result 



— ^ e H — -e 



aR 



R,) 

aR 



^3 = ^ V 



(A.3) 



aa aa 

with (T*j the spin connection for the constant curvature transverse space, satisfying d9^ + j A 
6^ = 0. The nonvanishing components of curvature form = (ia;°^ + u"'^ A cu^^ are given by 
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a a 



Rt 



a 
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e°Ae\ 
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e*^ A + 



J_ fR^ 
aR \ a 



R ra 



aa^ R 



Ae* 



1 [Rr 



RfO r 



a R \ a J ^ aR 



1 (R,r 



aR \ a 
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A e\ 



Rl_ 
a R 



Rr 



where TZ^ • is the curvature form of the constant curvature space 



(A.4) 



(A.5) 
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The components of the Ricci tensor in the tetrad frame Rab can be read from this expressions 
using 



Rab = acb. where n\ = A e'. 



(A.6) 



Then the Einstein tensor in the orthonormal frame Gab = Rab — \VabR, with R = ri°-^Rabi is 
given by 



G 



00 



Gil 



01 



{D-2 

{D-2){D-3) 



_^fRA ^ atR,t 
aR \ a / „ a^aR 



+ 



1 fRt 



R,t \ ^ / R.r \ ^ 

aR 



RrCl , 



aR) ^ i?2 



aR \ a J ^ aa?R 



(D-2)(D-3) 



aR) [aR) 



(A.7) 



{D-2) 



aR \ a J ^ aa^R 



1 / R^r \ ^ (^,rR,r 

aR \ a J J. aa^R 



(D-S) 



1 fR 



+ 



1 1 far 
— (a ta) ^ H ( — 



aa 



aa \ a 



P-3)p-4) 



0',tR,i 



aR \ a ) ^ a'^aR 



[aRj [aRj 



The next step is to compute the right-hand side of the Einstein equations. Again we work 
in the tetrad frame where the four-velocity is given by 



(A.8) 



so the energy- momentum tensor has the simple form Tab = diag {p,p, . . . ,p). With this we can 
write the full set of Einstein equations Gab = SnGNTab- For the sake of simplicity, however, it 
is convenient to replace the more complicated ij components by the Bianchi identities (12.51) . 
For the 00, 01 and 11 component of the Einstein equations one arrives after simple algebraic 
manipulations to Eqs. (12.41) . 
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Appendix B. Calculation of the coefficients in Eq. (12.271) 



In this Appendix we outline the procedure to compute the second coefficients C^^^(?/sp) of the 
expansion of (log?])' as a series expansion in powers of log{z/zsp). Taking a derivative with 
respect to log z in the last equation in fl2.22p and defining the right-hand side of the expression 
using the I'Hopital rule we arrive at 



2(v;^4(iogr;); = ;c;-(iogr^y(v;2)", 



(B.l) 



where by /C we denote the numerator of the right-hand side of the last equation in (12.220 



/C = 



d-2 

At the sonic point this satisfies 



■r]k+iS 



4-2d 
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As we have already argued, all derivatives at the sonic point should be a function of just 
one parameter that we can take to be ?/sp- Indeed, all the derivatives appearing in Eq. (1B.3P 
can be written in terms of y^p and (logr7)gp by using 
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together with 
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and 
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In addition, Eq. fl2.50p with the choice of the minus branch is also needed. 

Plugging all these expressions into Eq. ( IB.ip we arrive at a linear equation for (log?7)sp that 
determines also cip{ysp). The fact that this equation is linear is of conceptual importance. It 
means that, fixed the ambiguity in the choice of the branch determining [logr])'^^ there is no 
further ambiguity in the calculation of (log?7)gp. In fact, this extends to the calculation of all 
higher derivatives. The algebraic equations giving the values of the n-th derivative at the sonic 
point are always linear. 

Appendix C. Some details of the calculation of per- 
turbations 

In this Appendix we give the expressions for the derivatives of the coefficients needed to solve 
the perturbation equations around the singular points. In order to evaluate (13.161) the following 
derivatives are required 
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In addition, the derivatives of V^, logy and log// are also needed. For the first two we have 
whereas (log?])' can be written as 



(log?7)' = -T-v/r^, (C.3) 
where T and S are given respectively by 



+ 4rf - 6 - 2A;(d -l) + {d-2){k- 5)y, 

k — 1 

'^^tllls = 2(A; + 1)2(2/})^ f-V^' (C.4) 
y-1 \r] J 

- {d-2) \{d - 1)(A; - 2) + (3d - 5)y 
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